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Summary: The paper discusses the influence of wire length on the characteristics 
of a turbulent flow as measured with a hot-wire anemometer. Most of the mathematical 
treatment can be directly applied to other problems of length correction, such as may 
occur, for instance, in problems of astrophysical turbulence, as well as to some problems 
related to measurements of random processes. The discussion will, however, use the 
example of a hot-wire anemometer with particular attention to the influence of the 
length of the wire on the measurement of the intensity of turbulence, correlation 
coefficients, and scale and microscale of turbulence. 


Notation 


The use of a bar over a symbol denotes a time average taken over a time 
interval of at least 2T as defined in equations (8) and (9). 


Subscript “b” is used to denote the values determined by means of a wire of 
non-negligible length. 


Subscripts A and B refer to the wires A and B (see Fig. 3). 


The prime notation denotes a turbulent component and is also used in Section 4 
to designate the derivatives (such as R’, R”). 


O,X,.Y,Z rectangular co-ordinate axes (see Fig. 2) 
V(X, Y,Z,t) instantaneous vectorial velocity at time ¢ 
V’(X,Y,Z,t) turbulent vectorial velocity at time ¢ 
O,x,y.z rectangular co-ordinate axes with Ox parallel to V (see Fig. 2) 


u,v,w instantaneous velocity components in the directions Ox, 
Oy, Oz 


u’,v’,w’ turbulent velocity components in the directions Ox, Oy, Oz 
d,i,R, Ra, Ro, 9 
6,,A,B,V,a2,A, ~ quantities defined where used in Sections | and 2 
B,, E,v,¢,®,x,T 


*Some of the results given here were included in the United States Naval Ordnance Laboratory 
Memorandum 9658 (1948) distributed to a limited number of readers. 
The preparation of this paper was supported by the Bureau of Ordnance, Department of 
the Navy, under Contract NOrd 7386. 
Originally received April 1952. 
{The Aeronautical Quarterly, Vol. V, May 1954] 
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fluctuating voltage at the output of a hot wire amplifier 
constant defined by e= Kw’ 
length of the wire 


k constant defined by K=kl 
variables of integration 


R.(r) longitudinal correlation coefficient between the simultaneous 
turbulent velocity components u’ at two points on a line 
parallel to Ox and a distance r apart 


R,(r) transverse correlation coefficient between the simultaneous 
turbulent velocity components u’ at two points on a line 
parallel to Oy and a distance r apart 


L,,L, longitudinal and transverse scales of turbulence 


I, ratio of the correct to the measured longitudinal intensities of 


turbulence (/,?=e?/e,? is the ratio of the correct longitudinal 
turbulent energy to the energy measured with a wire of 


length /) 
A microscale of turbulence 
G factor defined by equation (35) 
4,5,P,Q defined with reference to Fig. 3 


Ry,» longitudinal and transverse correlation coefficients as 
measured with hot wires of length / 


Lyx, Ly, longitudinal and transverse scales of turbulence as measured 
with hot wires of length / 
Ky..Ky, ratios L,/L,. and L,/Ly,y respectively 


A», microscales of turbulence found from measurements, with hot 
wires of length /, of the correlation coefficients Ry and Ry, 
respectively 


1. Generalities of the Method of Measurement 


Let a fine metal wire of diameter d and length /, heated by an electrical current 
of intensity i, be placed normal to the direction of the velocity V of a gas flow 
(Fig. 1). If the velocity of the gas is constant, a state of equilibrium is reached 
between the heat supplied by the electrical current to the wire and the heat lost by 
convection, radiation, and conduction. The difference between the temperature of 
the wire 6 and the temperature of the ambient air 6, then remains constant. 


The relation between these quantities is given by the equation of L. V. King"’ 


?R=(A+B/V)(6—4,), 


(1) 


where A and B are factors depending on the flow conditions, the size and the 
material of the wire, and the nature of the gas, and where R is the electrical 
resistance of the wire at the temperature 9. By expressing the temperature 


2 


il 
yf 


iS 


TURBULENCE INVESTIGATIONS WITH A HOT-WIRE ANEMOMETER 


Fig. .1. 
A platinum wire fixed at the end of two prongs. 
The size of the wire can be estimated by com- 
parison with the millimeter scale. 


difference as a function of the electrical characteristics of the wire, it is found that 
ai?R.R 
R-R, 
with R, and R, the resistances at 0°C. and at the ambient temperature, and 2 the 
temperature coefficient of resistance. 
Equation (2) gives a relation between the wind velocity, the intensity of the 
heating current, and the resistance of the wire. One of the last two parameters can 
be replaced by the potential drop at the terminals of the wire, E=iR. 


If now the resistance R is kept constant, then the potential drop is given as a 
function of the velocity by the relation 


=A+B/V ; (2) 


R(R—R,) 


When, on the contrary, the intensity i is kept constant, the voltage is given by 
(A+BVV)R, 


The relationship between the potential drop and the velocity might be 
modified by including ‘in the electrical circuit a non-linear amplifier. It is possible 
to obtain in this way the relation 


where A, and B, are factors depending on the equipment and the flow conditions. 
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In each of these three cases the voltage is given as a function of the velocity 
and may be used to measure it. In general the function 


will be found experimentally by slowly varying the hamenend and measuring the 
voltage. 


When the direction of the velocity makes an angle # with the wire, and if this 
angle is not much different from «/2, the value of the potential drop is given as a 
function of the component of the velocity normal to the wire, namely 


For the experiments of G. B. Schubauer and P. S. Klebanoff?’ equation (7) is true 
for values of ¢ as small as 20°. 


When the velocity of the flow changes, its fluctuations are not followed 
instantaneously by the fluctuations of the voltage. There exists a delay due to the 
thermal lag of the wire. A so-called time constant of the wire describes the 
magnitude of the time delay. This time constant decreases with the diameter of 
the wire®’. The influence of the thermal lag becomes important when the wire is 
used in a flow for which fluctuations of velocity are of high frequencies. The 
influence of the lag can be compensated by placing a special electronic stage in the 
electrical circuit of the wire. The role of the compensation becomes essential when 
the hot wire is used for measurements of the characteristics of turbulence. The 
present paper does not deal with the question of compensation, and it is assumed 
here that the hot-wire equipment has no thermal lag or that its lag is perfectly 
compensated by an electronic stage. 


While we consider here the wire by itself, it should be remembered that it is 
connected to a heating source. The fluctuations of the velocity changing the 
electrical conditions of the wire have, therefore, an influence on the heating circuit. 
For instance, if it is desired to operate at constant current, we allow the resistance 
of the wire to fluctuate with the velocity of the gas. The overall resistance of the 
heating circuit of which the resistance of the wire is a part will be affected by these 
fluctuations”). However large the voltage of the source and however small the 
resistance of the heating current, there may still remain non-negligible fluctuations 
of the intensity. In the present paper this effect is neglected. 


2. Measurement of Turbulent Velocities with a 
Hot-wire Anemometer 


Let V(X, Y,Z,t) be the instantaneous vectorial velocity of a turbulent fluid 
at the time ¢ and at a point P(X, Y, Z) defined with respect to stationary orthogonal 
axes OXYZ (Fig. 2). Let us define the mean velocity at this point by 


t+T 


+ 
V(X,Y, 2)= lim [ Y,Z,s)ds, . ‘ (8) 
T 


2T 
t— 
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where the averaging is made over the time. We assume that, at any time f, there 
exists at every point of the turbulent field a sub-interval of time 2T such that for 
any interval of time equal to or larger than 2T the difference 


t4T 


t-—T 


is negligible. We assume, moreover, that the experimental conditions are such 
that we can consider the gas as being a continuum with a defined instantaneous 
velocity at every point. These assumptions are acceptable for most of the cases 
where a hot-wire anemometer is used, for instance, in ordinary wind tunnels. 


The instantaneous vectorial velocity at the point P(X,Y,Z) can be divided 
under these conditions in two parts: V(X, Y,Z) the mean vectorial velocity and 
V’ (X,Y, Z, t) the turbulent vectorial velocity, and therefore 


Let there be a second set of orthogonal axes Oxyz with the axis Ox in the direction 
of the mean velocity V(X.Y,Z). If u, v, w are the components of V (x, y, z, t) 
parallel to the axes Ox, Oy, Oz and if for each of these components their mean and 
their turbulent parts are found, then 


u(x, y,z.=u(x, y,z)+u' (xy, Z. 0) 
v (x, y, Z, (x,y,z. 0) 
w (x, y, Z, w’ (x, y, z, 0) 


v =w=0 due to the choice of the direction of Ox and where 
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There will be around the wire a spatial domain in which the wire integrates, 
during the measurement, the effect of the velocities at each point of the domain’. 
A rough approximation would be to consider that this domain is equal to the 
volume of the wire and that the weight factor of each point at this domain is the 
same. Whatever is the real spatial domain and the distribution of the weight factor, 
it is assumed in this paper that the turbulence, as well as the mean velocity, is 
homogeneous in this spatial domain. If, for instance, ® is an instantaneous value 
of any fluctuating physical quantity, we would have throughout the domain 


concerned 


The characteristics of turbulence do not change therefore by a translation 
of the axes. 

Consider now a wire placed in a turbulent flow normal to the mean velocity. 
In the co-ordinate system, the axis Ox is parallel to the direction of the mean 
velocity and the axis Oy parallel to the axis of the wire. It is assumed in the present 
section that the dimensions of the wire are negligible and that it measures the 
velocities at a point. It is assumed also (through the whole present note) that the 
turbulent velocities are small in comparison with the mean velocity, i.e. 


When the angle ¢ is approximately «/2, equation (7) can be applied. The instantan- 
eous voltage will be, therefore, equal to 


and if we neglect w”” in comparison with u? it becomes 


After expansion in a series and neglecting all the terms except the first two, the 
last equation gives for the fluctuating voltage 


If both sides of (14) are averaged, the mean voltage obtained is 


Under the assumptions of the preceding section, /(u) is independent of the time. 
The turbulent voltage is found by subtracting the last equation from equation (14), 
,av{u) 


For the measurements of turbulence, the electrical circuit includes an amplifier 
which, besides compensating the thermal lag of the wire, eliminates the mean 
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voltage E and amplifies the turbulent voltage E’. An ideal amplifier for the hot- 
wire anemometer, adapted to the explored turbulent flow would: 


(a) compensate exactly the influence of the thermal lag of the wire, 


(b) eliminate the mean voltage taken over a time interval at least equal to 
the time interval over which the average value for the fluctuations of the 
velocity should be computed, and 


(c) have a constant factor of amplification and no phase distortion for the 
entire range of frequencies of the turbulent fluctuations. 


At the output of such an amplifier we would measure a fluctuating voltage 


: . (18) 


where x would be the factor of amplification. 


e=xE’= x 


The derivative dy (u)/du is, in general, a function of the mean velocity u. As 
it was assumed that every experiment is made in a field of constant mean velocity, 
this derivative is considered as remaining constant during the measurement; its 
value will, however, change when u is changed. Equation (18) can be written 


where K = xdv) (u)/du is constant for each measurement. The object of the present 
paper is to discuss the influence of the length of the wire on the measurement of 
the characteristics of turbulence. In the following sections, after assuming that 
equation (19) would be verified for a hot-wire of negligible length, a corresponding 
relation will be written for a wire of non-negligible length. The equations will then 
be found representing the influence of the length of a hot wire on the measurement 
of various characteristics of turbulence. 


3. Measurement of the Longitudinal Energy of Turbulence with a 
Hot Wire Normal to the Direction of the Mean Velocity 


Let there be a fine wire of negligible diameter and of length / placed in a 
gaseous flow normal to the direction of the mean velocity. The turbulence of the 
flow is supposed to be homogeneous and the mean velocity constant in the domain 
of space-time in which the measurement is made. It is assumed that the thermal 
lag of the wire is perfectly compensated and that the amplifier eliminates the mean 
voltage and amplifies by a constant factor the fluctuating voltage for all frequencies 
included in the turbulence. Consider the orthogonal axes Oxyz with Ox parallel to 
the direction of the mean velocity and Oy superposed on the wire. The fluid is 
assumed to be a continuum with defined instantaneous velocities at every point, 
continuously differentiable with respect to each of the variables x, y, z and ¢. If 
under these conditions, an element of length dy of the wire is considered, then, at 
the time ¢, it is responsible in the overall output voltage of the amplifier for a part 
given by 

de 


dy Ku . » (20) 
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with K constant during the experiment. It should be pointed out that equation (20) 
is, in fact, an approximate relation and its use for a given experimental condition 
has to be justified. The factor K should be determined from experimental 
calibration curves. 


If the length of the wire is very small in comparison with the scale of turbulence 
and if the variation of uv’ along the wire is neglected, then equation (19) is obtained 
and is rewritten as 


where kl=K is used to simplify some of the equations which will follow. In the 
last equation, as in those which follow, the influence of the terminals of the wire, 
where it is fixed to the supporting prongs, is neglected. 


When the length of the wire is not negligible (compared with the scale of 
turbulence), the variation of u’(t) along the wire cannot be ignored. If e, () 
represents the instantaneous fluctuating voltage measured with such a wire, then 
by integrating (20) we find 


° 


A thermic ammeter measures fluctuating voltages in terms of root mean square 
values. The time interval over which the averaging is made depends on the 
characteristics of the ammeter. It is assumed that this time interval is at least equal 
to the time interval for which a constant mean velocity may be defined. Various 
other integrating apparatus beside the thermic ammeter might be used and a 
regulable averaging time might be provided. 


If such averaging apparatus is connected through the amplifier to a wire of 
negligible length, then one measures 


and for a wire of non-negligible length 


After transforming the integral given in the last equation into a double integral and 
remarking that the average of this integral is equal to the integral of the average, 


it is found that 
l 


l 


In the last equation there will now be introduced the transverse correlation coeffi- 
cient which is defined in a homogenous turbulence by the relation 


u’ (x, 0, z, t) u’ (x, y, Z, 


R, (y)= (26) 
u 
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with u?—[w (x,0,z. =[w’ (wy, z. OF. 


Equation (25) in which the coefficient of correlation is introduced, is now 
written 


If now 
and taking into account that R, is an even function, then 
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9 
For each of these double integrals the order of integrations is inverted. After 
finding the new limits of integrations, one obtains 

l I l-s 


ds R, (s) ds, + ds R, (s) ds, | 


° ° 


which leads to 2k*u? | (/—s) R, (s) ds. (28) 


The last equation gives the value of e,? measured with a wire of non-negligible 


length /, as a function of the true value of u’* and of the transverse correlation curve 
R,(y). Comparing this equation to (23), and defining by /,? the factor by which 
the measured longitudinal turbulent energy* e,? should be multiplied to give the 
correct value of this energy, it is found that 


where 1,2 =e" / . G0) 


As R,(y)<1, the turbulent energy measured with a hot wire of non-negligible 
length is smaller than the true value of this energy. Equation (29) was first found 
by H. K. Skramstad® by a different method and the present author showed it by 
still another method". 


3.1. Isotropic turbulence 

If the turbulence is homogenous and isotropic, then instead of expressing /, as 
a fraction of the transverse correlation coefficient R,, it is possible to express it as 
a function of the longitudinal correlation coefficient 


u’ (0, y, ) y,Z,0) 


R, (x)= (31) 


*We call longitudinal energy of turbulence the quantity 4pu*, where p is the density of the 
gas. By analogy we may speak about e,* as the measured energy and e* as the real energy, 
because the values of these measurements lead to a comput ipu u?, 
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The relation between the two correlations is given by the equation of Th. von 
Karman“? 


R,()=R:; 


(32) 
where r is the distance between the two points at which the turbulent velocities are 
considered. The distance r represents here x or y, according to the case. Eliminat- 
ing R, between (29) and the last equation, we find the factor of length correction 
I, expressed as a function of the longitudinal correlation by 


7) Reds. 


Integrating equation (32) from zero to infinity with the condition 


lim rR, (r)=0, 
1-00 


gives L,=2L,, 


co 
where i= | R..(s) ds, and 


are the longitudinal and transverse scales of turbulence. 


3.2. Approximate relations for large |/Lw 
(i) Considering the case when //L, is large (and neglecting yR, (y) and R, (y) 
for y</), one finds with equation (28) that 


where L,,,, is the first moment of the transverse correlation curve 


oo 
Lyy= { SR, (s) ds. 


° 


It appears, therefore, that when //L, is large, the longitudinal energy of the turbu- 
lence measured in the flow is a linear function of the length of the wire and decreases 
with increasing /. The factor J, can be given for this case by the relation 


where L,,,,/L,? depends on the shape of the correlation curve. When //L, increases 
indefinitely, it is found at the limit 


lim (//1,?)=0, ‘ (38) 
1/L 00 
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and lim (e,?/k7/?)=0. 
1/L 300 


The longitudinal turbulent energy measured by a wire whose length is very large 
in comparison with the transverse scale of turbulence will approach zero even if 
the real energy is not negligible. 


(ii) In homogenous and isotropic turbulence we can apply equations (33) and 
(35). For large 1/L,, when //L,=1/(2L,) is also large, we find 


= IL, =2k?u"? IL,, 
(40) 


At large values of //L,, the energy of longitudinal turbulence, measured with a wire 
of non-negligible length in a flow of isotropic turbulence, is inversely proportional 
to the length of the wire. 


3.3. Approximate relations for small 1/L. 


(i) Take now the case when the length of the wire / is small, compared to 
the transverse scale of turbulence L,. One will have to consider now the correlation 
coefficient R,(y) for small values of y only. Developing in a Taylor’s series, 
we find 


= 


If all terms of the series except the first three are neglected and if R, (y) in equation 
(29) is replaced by its value, then . 
1 P @R,O)  dR,@) 


(ii) Consider now the case of isotropic turbulence with the two correlation 
coefficients R,(r) and R,(r) developed in series similar to (41) and introduced in 
equation (32). Neglecting again the higher order terms we shall find 


(0) 


(0) 
dr? 


R,()=1+7 


+47 (43) 
Using the definition given by G. I. Taylor’ for the dimension A of the smallest 


eddies responsible for dissipation of the turbulent energy 


1 d?R,@) (44) 


and the notation of G. K. Batchelor and A. A. Townsend"° 
d‘*R, (0) 
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one finds for the transverse correlation coefficient the relation 
G 
R,(y)=1- +$ 
Using the same notations, equation (42) will give 


Me 
120 


The last equation could also be found directly from (33). 


In (47) the factor J, is given as a function of //A. This equation was found 
for small values of //L,. It is seen now that its use should be considered only for 
the cases when //A<1. For very small //A it might be possible to use even the 
more simplified relation 


1 P | 


It should be noticed that in both equations (47) and (48) the factor J, is given as a 
function of the true dimension A and not A, as it would be measured by a hot wire 
of non-negligible length. The relation between the true value of A and the value of 
A, found from measurements with wires of non-negligible length will be considered 
later. The coefficient G depends, as A does, on the shape of the correlation curves 
at small values of r. 


3.4. Measurement with a wire whose length | is of the 
order of magnitude of Lv 


When the length of the wire cannot be considered as very large or very small 
compared with L,, then it is necessary to use the general equation (29). In that 
case the factor J, depends on the shape of the correlation curve. 

(i) Assuming a transverse correlation function to be of the form 


which was used also by Skramstad’, it is found that 


[ ( | 


(ii) For a transverse correlation function 


R,0)=exp(~ 73), 
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the value of J, is given by 


erf (a)= exp (—s”) ds. 


(iii) If now the turbulence is isotropic and the longitudinal correlation function 


is given by 


then, using equation (33), one finds 


ert (4 i): 


or if equation (35) is applied, the factor /, might be expressed as a function of the 
transverse scale of turbulence by 


2L, Jr 


4. Measurement of Correlation Coefficients with Two Hot Wires 
Normal to the Direction of the Mean Velocity . 


Two identical wires A and B each having length / and a negligible diameter are 
placed in a plane making an angle 4 with the direction of the mean velocity U 
(Fig. 3).. The two wires are perpendicular to the direction of the mean velocity 
and are placed between two planes which are parallel to the direction of 
U and perpendicular to the plane A. Let r be the distance between the two 
wires, x the projection of r on the direction of the mean velocity, and y its projec- 
tion on a plane perpendicular to this direction. The two wires are heated to the 
Same mean temperature and the electrical equipment is built in such a manner that 
to the same fluctuations of the velocity correspond for each wire the same voltage 
fluctuations. A switch allows the instantaneous potential drops at the terminals of 
the two wires to be added or subtracted. The resulting voltage is fed into an 
amplifier which eliminates the mean voltage and amplifies the turbulent voltage. 
The thermal lag of the wire is assumed to be perfectly compensated. The output 
of the amplifier is then sent to a meter measuring root-mean-square values of the 
fluctuating voltage. 


Consider first the case when the length / is negligible compared to the transverse 
scale of turbulence. If the instantaneous fluctuating voltages, e,, es, at the terminals 
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Fig. 3. 


of the two wires, A and B, are added, then the averaging meter measures the value of 


where u,’ and u;’ are the longitudinal turbulent velocities at the two wires. If the 
voltages, e4, es, are opposed, the measure will give 


— és) (u4’ — Uz’). ° (57) 


Dividing the difference between these two measures by their sum, we obtain 


In a flow of homogeneous turbulence 

and, in consequence, 


Let us represent by R, (r) the coefficient of correlation along the direction A (corre- 
lation between the longitudinal turbulent velocities at two points placed at a 
distance r apart on a straight line which makes an angle A with the direction of the 
mean velocity). This coefficient is found’, with equations (58) and (60), as 
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Consider now the case where the length of the wires is not negligible in com- 
parison with the scale of turbulence. If R,, , (r) is the correlation coefficient along 
the direction 4 measured in that case and, ¢),,, év,2, are the instantaneous potential 
drops at the terminals of the wires, then 


+ 


(62) 


R, = 


Let P and Q be two points placed on the wires A and B, respectively, at the distances 
Ss, and sz; from the chosen plane of origins (Fig. 3). The instantaneous potential 
drops, éy,,4 and é,, at the terminals of the two wires are given by the relations 


u’ (s4, ds, and u’ (Sz, t) dsx. (63) 


Since the turbulence is homogeneous, 


[u’ (s4)}? =[u’ 
and = = 
Equation (62) can then be rewritten 


Cv,ACb,B 
b R, r)= 


Applying equation (63), we find 


H a8 


=k? (S4,t) U’ (Sg, t) dsp =k? | (S4,t) (Sp, t)dsadsy.. (65) 


Let R; be the correlation coefficient between the longitudinal turbulent velocities 
at the two points P(s,) and _— The distance PQ is equal to vy {r?+(sg—s,)*} 
and we have 


(S4,t) (Sp. t) 
u? 4 


(66) 


Rs( = 


which will give with equation (65) the relation 


2 


The function Rs is a scalar function of a scalar and a vector. To perform the 
integration, this function has to be transformed into an ordinary scalar function. 
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4.1. Isotropic turbulence 


The transformation of Rs into an ordinary scalar function can be easily made 
when the turbulence is homogeneous and isotropic. Applying the equation of Th. 
von Karman, we can express R; as a function of the longitudinal and of the 
transverse correlation coefficients R, and R, by 


= Re (vir + (Sx — 84)? 


y? 
r’? 


Ry (V (68) 


If Rs is replaced in equation (67) by this value, then the following expression 
is found 


l 1-8, l 


ku | ds, R,.(V ds+ ds, (J |, 


(69) 


where S=Sp— 


Both coefficients R, and R, are even functions of /(r?+5”), and the integrals of 
the last equation can be transformed by operations similar to those which were 
used in Section 3. After changing the order of integration in the double integrals 
and finding the new limits of integration, we obtain 


l 


The correlation coefficient R, , measured with the wires of non-negligible length 
can now be obtained as a function of the real longitudinal and transverse correlation 
coefficients R, and R, from equations (23), (64), and the last equation. We find 


(70) 
where the factor /, is given by equation (30). 


(i) It is, of course, possible to eliminate one of the functions R, or R, from 
equation (70) by applying the von Karman equation (32). 


(ii) When the two wires are placed in a plane parallel to the direction of the 
mean velocity (A=0°) with their axes perpendicular to this direction, then y=0 in 
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equation (70). Eliminating now R, between equations (32) and (34), we express 
the measured longitudinal correlation as a function of the real longitudinal 
correlation function by 


where J, is given by equation (33). 


4.2. _Axi-symmetric turbulence 


The case of axi-symmetric turbulence with the axis of symmetry parallel to 
the direction of the mean velocity includes the special case of isotropic turbulence. 
The relation found hereafter is, therefore, also valid for isotropic turbulence. For 
the measurement of the correlation coefficient, the two wires are placed in a plane 
perpendicular to the direction of the mean velocity (A=/2). In this case x=0 in 
equation (70), and we find for the measured transverse correlation coefficient 


Roy (y)= | Ry (v ds, (72) 


with the factor J, given by equation (29). This equation was found first by H. K. 
Skramstad“ by quite a different method. 


In equations (70) to (72) the correlation coefficients, measured with hot wires 
of non-negligible length / are expressed as functions of this length, of the correct 
correlation function, and of the length correction factor /,. The factor /, is itself 
given as a function of / and is a function of R(r) through equations (29) or (33), 
according to the case. In order to find R,(r), it is necessary, therefore, to know 
the form of the function R(r). The shape of the correlation curve varies with 
the experimental conditions, and various equations will have to be used to 
represent this curve. It is often-difficult to represent the correlation function by 
an equation which gives values of R(r) in good agreement with the correlation 
coefficients over a large interval 0, r and is simple enough for application in 
equations (70) to (72). It is in general easier to represent the experimental correla- 
tion curve by two equations: one for small values of r and another for r larger 
than some given value r,. The factor /, is given as a function of R (r) in the interval 
0, 1, while the integrals of equations (70) to (72) are functions of R(r) in the 
interval r, /(r?+/?). We should try, therefore, to use an equation which accurately 
represents R (r) in the interval 0, / to compute /, and an equation accurately repre- 
senting the correlation curve in the interval r, /(r?+/*) for equations (70) to (72). 
It will evidently be preferable to represent the correlation curve by a single equation 
if it is possible to do so without excessive difficulties. 
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4.3. Approximate relations for large r/l 


When the length / of the wire is small compared with the distance r between 
the two points at which the longitudinal turbulent velocities are considered and 
this distance is not small compared with the scale of turbulence, then equations (71) 
and (72) can be simplified. As the interval of integration is 0<SsS l, we find by 
neglecting the terms with s*/r* and higher that 


V(r? +s?)= 


Developing R (/ {r?+5?}) in a Taylor’s series and neglecting the terms with s*/r* 
and higher, we find 


For the case of isotropic turbulence, equations (71) and (73) give the measured 
longitudinal correlation coefficient 


where the factor /, is given by equation (48) or (49) according to the magnitude 
of //A. The last equation can also be written as follows 


In this equation the correct coefficient of correlation R, is given as a function of 
R,»,. measured with hot wires of non-negligible length and of the values of R,’. As 
a first approximation, we can use in equation (75) the measured values R,,.’ (x) 
and R,.” (x) instead of the correct values R’ (x) and R” (x). 


For very small values of //x, the terms with /* can be neglected. At very large 
distances x, when //x is still smaller, we shall have simply 


In a simpler way for the case of lilaiaatie turbulence equation (72) gives 
with (73) for the measured transverse correlation coefficient 


(77) 
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The last equation might be used to correct the measured values of the transverse 
correlation coefficient if, as a first approximation, R,,,’ (y) and R,,,” (y) are used in 
(78) instead of R,’(y) and R,” (y). 


When //y is sufficiently small we can use the simplified equation 
1 
Ry 73 Row 0). 


4.4. Influence on the shape of the measured correlation curve 


The value of Ry, as measured with hot wires of non-negligible length, can be 
computed if we assume some given function to represent the correct correlation 
curve. 


In axi-symmetric turbulence with a transverse correlation function of the form 


R, (y)=exp(—| y |/L,). + (80) 


the corresponding R,,, (y) values were computed numerically by H. K. Skramstad“”. 
When the correlation function is of the form 


we find by applying (72) 


Replacitig J, by its value given by equation (52), gives 


Consider now the case of isotropic turbulence with the longitudinal correlation 
function 


tx 
R.(x)=exp(- x). 
Eliminating R, and J, between the last equation and equations (35) and (54) gives 
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The transverse correlation function corresponding to (83) is found with G2). 


and applying equations (72) and (55), we find for the transverse correlation 
measured with wires of non-negligible length 


1- exp(- 
Bo 


5. Measurement of the Scale of Turbulence 


The scale of turbulence (longitudinal or transverse) is found by taking the 
integral of the corresponding correlation curve from zero to infinity 


oo 


= | R (s) ds. ‘ (87) 


In general, the value of L cannot be obtained with great precision because the 
upper limit of integration is infinity. 


When the length of the wires is not negligible, the measured scale of turbulence 
is equal to 
ce 


The factors by which the measured values of the longitudinal and transverse scales 
of turbulence L,, and L,,, should be multiplied to give the real scales of turbulence 
L, and L, are 


(i) For axisymmetric turbulence with a transverse correlation function of 
the form 


Ky = (89) 


R, (y)=exp(- | | /Ly), 
it is found” that 


Kyy=2h? )]. 


where J(a~)= + exp (—acosec dB, 


and where J, is given by equation (50). 
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(ii) When the transverse correlation function is 


R,Q)=exp(- 73). 


we find that Ky,=1 and . 


(iii) If we consider the case of isotropic turbulence with a longitudinal 
correlation function 


R. (x) 


then K,,=1 
Applying equation (50) we find for this case that 


!-exp(- 
= 


6. Measurement of the Microscale of Turbulence 


The Taylor dimension A of the smallest eddies responsible for the dissipation 
of the turbulent energy, which is often called the microscale of turbulence, is a very 
important characteristic of isotropic turbulence. The dimension A is much better 
defined theoretically than the scale of turbulence, being given simply by the relation 


R.” (0) (0) 


(95) 
From equations (42) and (43) we find for small distances r the relations 

(96) 
(97) 


of which the second was given before as equation (46). We shall consider the two 
cases: (a) when the microscale of turbulence is obtained from the longitudinal 
correlation curve and (b) when it is found by using the transverse correlation curve. 
It should be remembered that in the present paper we consider only the influence 
of the non-negligible length of the wires on the measurements, assuming that from 
every other viewpoint the hot-wire equipment is perfect. For instance, we neglect 
for the present the effect of the wake of one wire on another, which may be quite 
important when R, for small x is measured. We neglect here also the effect of 
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radiation between the wires, and of the perturbation of the flow, which may not be 
negligible when the distance between the wires is very small. 


Consider now the case when the microscale of turbulence is measured with 
two hot wires of non-negligible length placed in a plane parallel to the direction 
of the mean velocity. The measured A,,, will be defined by the relation 


1 
Anz R,,.” (0) ° (98) 
From equation (71), 
l 
1 


If the microscale of turbulence is measured with two wires of non-negligible 
length placed in a plane perpendicular to the direction of the mean velocity, then 
the measured A,,, is defined by the relation 


2 
Any Mig (0) > (100) 
and, from equation (72), 
l 
1 I,” (I 5) 


When the turbulence is assumed to be isotropic, then equation (32) can be 
applied and we obtain, with equation (101), 


In this relation the microscale of turbulence measured with two parallel wires placed 
in a transverse plane is given as a function of the longitudinal correlation 
coefficient R, (x). 


We find, moreover, for the case of isotropic turbulence the relation 


Aas” 1 
— 


(103) 


> 
wo 
N| 


R. (5) ds 


= 


where the expression in the brackets tends to zero with /. The last equation 
shows that the value of the microscale of isotropic turbulence, measured with two 
wires of non-negligible length placed in a plane parallel to the direction of the 
mean velocity, is different from the microscale measured with wires placed in a 
plane perpendicular to the direction of the mean velocity. 
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When the length of the wire is small compared with A, then equations (96) and 
(97) may be used to represent the correlation coefficient. In what follows we shall 
neglect the terms with (//A)*. 


When A is obtained from measurements with two wires of non-negligible 
length placed in a plane parallel to the direction of the mean velocity, then we can 
apply equation (99). After finding R.’ from equation (96) and replacing 7, by its 
value given in (48), we have 


This relation gives the factor by which the measured value of A,” should be 
multiplied to give the correct value A*. We shall notice that A,.=A for G=3 
and that the real A is smaller than the measured A, for G > 3, which in practice 
seems to be the most general case. Instead of being expressed as a function of (//A) 
the correction factor (A/A,,.) can be given as a function of (//A,,.). Solving equation 
(104) for A gives 


If now the microscale of turbulence is obtained from measurements with two 
parallel wires of non-negligible length placed in a plane perpendicular to the 
direction of the mean velocity, then we apply one of the two equations (101) or (102). 
It will be found that 


(104) 


The real A is smaller than the measured A,, as long as G>4. If we express the 
correction factor (A/A,,) as a function of (1/A»,), we find 


The ratio of the microscale of turbulence measured in a plane parallel to the mean 
velocity direction to the one measured in plane perpendicular to this direction 


(106) 


is given by 
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When the distance r is small, we find (for isotropic turbulence), using equations 
(96) and (97), the following simple relation between the longitudinal correlation 
coefficient R(x) and the transverse correlation coefficient R, (y) 


1—R,(r) _ 

i- R, (r ) 
The similar relation for the correlation coefficients measured with wires of non- 
negligible length will be 


1 Row (r) 2 Ap” 


This ratio is larger than 2. 
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On a Simple Method for Calculating Laminar 
Boundary Layers 
K. E. G. WIEGHARDT 


(Formerly Admiralty Research Laboratory, Teddington, now at 
Hamburg University) 


Summary: A simple one parametric method, due to A. Walz and based on the 
momentum and energy equations, for calculating approximately laminar boundary 
layers is extended to cover axi-symmetric flow as well as plane flow. The necessary 
computing work is reduced a little. 


Another known method which requires still less computing work is also extended 
for axi-symmetric flow and, with the amendment of a numerical constant, proves 
adequate for practical purposes. 


1. Introduction 


Since there are already several methods of calculating approximately laminar 
boundary layers for incompressible flow, it would seem necessary to justify the 
development of yet another method. The following method tries to obtain in a 
practical manner, results as correct as possible and with the smallest amount of 
computing work. Obviously a very accurate method is to represent the velocity 
profiles approximately by a two parametric class whereby the partial differential 
equation for boundary layers is replaced by two ordinary differential equations: 
the momentum and the energy equation”. But here the interpolation in the two 
parameters requires considerably more computing work than that with the usual 
one parametric methods. Therefore A. Walz’) altered this method back again into 
a one parametric method using the energy equation instead of the boundary condi- 
tion for the second derivative 0?u/dy? on the wall. By this means he obtained better 
results than by the usual one parametric methods based on the momentum equation 
only. So he proved that it is more important to satisfy the energy balance on the 
average over the whole layer rather than the usual boundary condition for the 
curvature of the velocity profile at the wall. As the latter condition is too sharp for 
a limited variety of profiles, a wrong profile is sometimes selected by it. The method 
of A. Walz can be altered slightly as follows. It is extended so as to calculate 
boundary layers on axi-symmetric bodies (without incidence) as well as in two- 
dimensional plane flow. Also, the computing work is reduced. 


Notation 


x the axial co-ordinate 
r radial co-ordinate 
r, the co-ordinate of the body of revolution 


Paper originally received January 1952. 
[The Aeronautical Quarterly, Vol. V, May 1954] 
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5 the distance measured along the surface of the body 


n the distance normal to the surface 


ncos a dx 
u,v the normal velocity components at any point in the layer 


U _ the velocity at the edge of the boundary layer 
U, _ the velocity of the undisturbed stream 

p _ the density of the incompressible fluid 

the viscosity 

v =p/p=the kinematic viscosity 


5 the thickness of the boundary layer 


“)" an=the displacement thickness 
U/r, 


a, = | dn=the momentum thickness 
= | an=the energy thickness 
D= J =the dimensionless dissipation in the layer 


H = 6,/6,=a parameter characterising the velocity profile 
H;, 8; / 8, 
R, = U8,/v=the local Reynolds number of the boundary layer 
R = U, R’/v=the Reynolds number of the main stream 
R’ a characteristic length 
X = 6,R,=U5,?/v=a useful computing quantity 
~=f (3) =the shear stress at the wall=the skin friction 
n=0 
é(u/U) 
(n/8,) 
_ 9,7 
U \én’? 


). = R,=the dimensionless skin friction 


> 
ll 


) =another characteristic profile parameter 
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2. The Calculation Method 


As the equations for axi-symmetric flow include those for plane flow and are 
only a little more complicated, the method may be derived for the boundary layer 
along a body of revolution, the co-ordinates of which are x (along the axis) and the 
radius r,(x). In the case of plane flow all the equations are to be simplified 7 
putting r, (x)=constant —> oo. 


The equation of continuity is 
0 0 
+ (rv) =0 (1) 


and the boundary layer equation 


ou dU u 
From the latter, two ordinary differential equations can be derived by multiplying 
by r or ur and integrating through the layer. These are the momentum equation 


df wf 
0 0 


and the energy equation 


0 

yer 

The physical meaning of the latter equation is that the energy dissipated by the 

viscosity equals the loss of kinetic energy inside the whole layer. 


Equations (3) and (4) can be written more shortly by introducing the quantities 
defined. Then equation (3) becomes 


1 dU 2 dr, 


dx { 
ds ds r, ds 


} 
This is the usual differential equation for X or 6,, the momentum thickness. Another 


equation for the profile parameter H can be derived by combining equations (3) and 
(4) as follows: 


dH 1 dU 1 


= H,,(H — 1) 32°” 
where f= “dH,,/dH and g= dH,,/dH (7) 
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TABLE I 
A B Cc D E F G 
2 | 0-83 64 0-73 2:2 79 0:26 
| 4 0°25 6 9 26 
2:05 | —1:07 0:58 72 0:79 2:5 8-7 0:35 0-439 
| 4 0:27 09 7 0°4 0-10 25 
21 | -1-03 8-1 0:86 2:9 9-5. 0:45 0-414 
| 5 0°30 1-0 7 0:4 09 0°10 24 
2:15 | —0-98 0-01 9-1 0:93 33 104 0:55 0-390 
6 0°34 1-1 9 0-5 1:0 0-10 23 
2:2 —0-92 +033 102 38 114 0:65 0-367 
7 0°39 1:2 9 0°5 11 0-11 22 
2:25 | -—0-85 0-72 114 43 125 0:76 0-345 
| 8 0°45 3 0-10 0°5 1:3 0°13 21 
2:3. | -0-77 1:17 1-21 48 138 0-89 0-324 
9 0°51 0-11 0°6 0°13 20 
2:35 | —0-68 168 141 1-32 54 153 1:02 0-304 
0-11 0°59 0°12 0:7 0°14 19 
24 |--0°57 2:27 15:7 1:44 6:1 169 1:16 0-285 
0°13 0°66 1°8 0°14 0:7 1:8 0-15 18 
2:45 | -0-44 2:93 175 1-58 68 18:7 1:31 0-267 
| 0-15 0°74 2:0 0°15 0:8 1:9 0°17 17 
2:5 | —0-29 3-67 19-5 1:73 76 206 1:48 0-250 
0°17 0°82 0°16 0-9 2°71. 0-18 16 
2:55 | —0-12 449 21-7 1-89 8:5 22:7 1:66 0-234 
0:18 0:90 2°4 0°18 1:0 2°4 0°19 15 
26 | +0-06 5:39 241 2-07 95 25-1 1:85 0-219 
0°19 0:99 | 0:20 1:0 2:6 0:22 15 
265 | 0:25 6:38 268 2:27 105 27:7 2:07 0-204 
| 0:21 1:09 2:9 0:22 0:23 14 
2:7 | 0:46 747 297 #249 116 306 2:30 0-190 
0-23 1:21 3:2 0:24 EZ 3:2 0:25 13 
2:75 | 0:69 8-68 32:9 2:73 128 338 2-55 0-177 
0:26 1°32 0:26 1:4 0:28 12 
28 0:95 10-0 364 2:99 142 37:4 (2:83 0-165 
0-29 3:9 0:28 40 0:29 11 
0:32 4:3 0°30 1:7 4:4 0-31 11 
29 | 1:56 13-2 446 3:57 174 458 3-43 0-143 
0°36 1:9 4:7 0°32 2:0 48 0°33 10 
2:95 | 1-92 15-1 49:3 3:89 194 506 3:76 0-133 
0°40 20 0°34 5°3 0°35 10 
| 2-32 17:1 545 423 216 559 411 0-123 


3-0 


These equations are still exact, of course, in general it would be impossible to 
compute the boundary layer only by equations (5) and (6), as these two ordinary 
differential equations cannot replace the partial differential equation (2) and all its 
boundary conditions. However, restricting oneself to an approximate calculation 
by assuming a one parametric class of velocity profiles, H,,, ¢ and D and conse- 
quently f and g become functions of the one parameter H alone: f(H) and g(H). 
In this case obviously equations (5) and (6) enable X (s) and H (s) to be computed 
if the body r,(s) and its pressure distribution, or U (s), are given. 
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For axi-symmetric flow this assumption of a one parametric class implies 
neglecting the term r/r, during the integration across the layer (see the definitions 
of the characteristic quantities). This is the same approximation Mangler “) made 
in his paper about the transformation of an axi-symmetric boundary layer into a 
corresponding plane one. Since usually the layer is very thin when compared with 
the radius of the body of revolution, the assumption r/r,=1 is good enough as r 
varies only between r, and r,+6. The influence of the shape of the body of revolu- 
tion still remains in the term (2/r,)(dr,/ds) X in equation (5), whereas equation (6) 
is approximately the same whether the flow is axi-symmetric or plane, as it is 
independent of r,. 


The system of the two differential equations (5) and (6) is to be solved step by 
step. For that the step As may be chosen so small that (1/U) dU /ds can be taken as 
constant through the integration step. Then the local pressure gradient is 
characterised by 


1 


Y= ~ U ds as (8) 
where the index 1 denotes the beginning and 2 the end of the interval. Practically 
this means that the steps As must be made as small as is necessary to draw the curve 
U (s) with reasonable accuracy. 


In the same way all the other functions are linearised over the interval. Hence 
from equation (6) 


x,= where AH=H,-H, 
AH + (f+4fOH)y 
of ; ag 
f= dH and 2 dH’ (9) 


Introducing this expression for X in equation (5) a quadratic equation for AH is 
obtained. Yet for small steps As, AH will also be small so that the term with 
(AHY can be neglected. Then 


A(1+0)+By+C 


y 
D+ Ey+F y+ 


ry, ds + 


(»=0 for plane flow), and A, B, ... G are the following functions of H: --. 
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A=g.B=(5 5 5 B,F= 


(11) 


In equation (10) the argument H in these functions is to be taken at the beginning 
of the step: A(H,), B(H,), ... G(H,). Once AH is computed AX follows from an 
equation derived from equation (5): 


_ +H, +3)y+20] X,/As 
1+4((4,4+4.4+ 3)y+ 20] As. 2) 


According to its definition X gives the momentum thickness 


8,2=Xv/U or vR-,/ 


where R=U,R’/v is the Reynolds number of the undisturbed stream and R’\a 
characteristic length. With 5, giving the size and H the shape of the velocity profile 
the boundary layer is known once X (s) and H (s) are calculated. 


To establish the functions A, B, ... G a certain class of velocity profiles has to 
be chosen. Walz” has shown that the results do not depend very much on the 
particular class if the energy equation is used together with that for the momentum. 
This independence obviously is fundamental for the usefulness of the whole method. 
So the Hartree-profiles which have been proved useful before may be taken again. 
Hartree“ tabulated the velocity profiles for the case U ~ s” where the boundary 
layer equation (2) can be tackled analytically. For each of his profiles (for various 
m) the integrations for 5,, 5, and 5, were carried out numerically and the functions 
A, B, ... G calculated (see Table I). It was not necessary to integrate the dissipation 
D to find g(H) after equation (7), as for these profiles 


(07u 
gQ@= ~Af with A= -— U (14) 


For, in this case, U ~ s", the velocity profile does not change its shape at all along 
the distance from the stagnation point s as long as m is constant. Hence dH/ds=0 
and, according to equation (6), g= -(1/U)(dU /ds) Xf. Further, as these are exact 
solutions the general boundary condition for equation (2) 


A= ~{ and by definition A= — = 


is satisfied, and so simply g= - Af. Yet it is just this boundary condition (15) which 
has been dropped in favour of the energy equation in the approximate method 
already described. Therefore, it will be fulfilled only approximately when the 
relation  ~ no longer holds. 
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Fig. 1. 
Howarth’s flow: U=1—s. Comparison of the dimensionless skin friction ¢ calculated by 
various methods. 
3. Examples 


To check the accuracy of the calculation method the plane boundary layer has 
been computed for the case U=1-s, which Howarth has treated analytically by 
series. The calculation begins at s=0-05 where 6,=0-162, X=0-162? (1—0-05)= 
0-0249 and H=2-77, «=0-172. The result as given in Fig. 1 (the dimensionless 
skin friction over the distance s) deviates very little from the exact solution. The 
separation point is only found by extrapolating as the range of Table I is not great 
enough. Namely, between H=3-480 and 4-031 (for the separation profile) no 
profiles are tabulated by Hartree so that an extension of Table I is impossible 
without computing more velocity profiles afresh. However, in practice the laminar 
boundary layer will turn into a turbulent one at some distance before separation. 
Hence it would not pay to strive after greater accuracy near to the separation 
point itself. 


For comparison Fig. 1 gives the separation point as found by various methods 
when starting at s=0-05: 

(a) Th. v. Karman-C. B. Millikan s,,,=0-102 

(b) A. Walz (Hartree profiles; based on the momentum equation only) =0-103 

(c) A. Walz (Hartree-profiles)®) =0-114 

(d) Two parametric class of velocity profiles) =0-116 

(e) H. Gértler (relaxation method) ) =0-118 

(f) L. Howarth (series)®)=0-120 

(g) D. Meksyn (series) =0-1235 

(h) Present method=0°125 

(i) A. Walz (profile: polynomial of the fourth degree)? =0-125 

(k) H. Schlichting (profile: polynomial of the sixth degree) =0-138 

() K, Pohlhausen (profile: polynomial of the fourth degree) =0-160, 
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Half-body of revolution (produced by a single source): Velocity outside the boundary layer, 
momentum thickness, skin friction and neutral stability point for various Reynolds numbers. 


Further, the boundary layer has been computed for three types of head of an 
infinite long body of revolution: 


(a) half-body 
(b) hemisphere and cylinder 
(c) }-calibre rounded head and cylinder. 


The pressure distribution in (b) and (c) was that calculated recently by Vandrey"*? 
by a new method (applicable to bodies with discontinuous curvature of the meridian). 
As these pressure distributions are calculated for ideal fluid the results as given in 
Figs. 2, 3 and 4 hold only for high Reynolds numbers. 


The strokes on the surface of the bodies indicate the position of the neutral 
stability point at the respective Reynolds number, i.e. the point after which the 
laminar layer might become unstable. The actual transition to the turbulent 
boundary layer will occur between this point and the separation point. The estima- 
tion of the neutral stability point is given in the Appendix. 


In (b), U (s) was given only at a few points s, as indicated by Fig. 3, so that the 
actual curve is rather uncertain near the end of the head. But obviously as the 
separation is there in any case it does not matter very much. 


Each case needed only a couple of hours’ computing work. 
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Hemispherical head: Velocity outside the boundary layer, momentum thickness, skin friction 
and neutral stability point for various Reynolds numbers. 


4. Comparison with Another Very Simple Method 


For plane boundary layers a somewhat simpler calculation method has been 
proposed by H. Holstein and T. Bohlen”*) and A. Walz"*) in Germany and by 
Young and Winterbottom"" and, in detail, by Thwaites"” in this country. It is 
based on the momentum equation (5) alone which might be written 


dX 


Supposing that the velocity profiles can be represented by a one parametric class, 
the right side of (Sa) is a function of the parameter H alone. For « this is obvious. 
Yet it holds also for (X /U) dU /ds which equals another parameter 


07u/U 

- 

(compare equation (15)) according to the general boundary condition which has 
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Fig. 4. 


1/4-calibre head: Velocity outside the boundary layer, momentum thickness, skin friction and 
neutral stability point for various Reynolds numbers. 


been dropped in the method already described. Hence, the left side of equation (5a). 
dX /ds, is also a function of H alone. 


Using Thwaites’s notation for A= — m, «=/(m) and L(m)=(2H + 4) m +21, the 
momentum equation becomes 


L(m)-—m. R . (16) 
ds 


For various known solutions of plane boundary layer problems the relation L(m) 
is almost the same, which proves the main assumption of a one parametric class of 
velocity profiles. Furthermore, L(m) can be approximated by a linear function 
such as 


(17) 


L(m)=a+bm. 


Then equation (16) can easily be integrated to give an explicit formula 
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and the parameter is to be found by 


_ XdU 
U 


(19) 


For axi-symmetric flow Mangler’s transformation? alters this equation (18) into 


x= | U'rtds 20) 
0 


with r, (s)=local radius of the body of revolution. This might also be written 


3/R’ 


m= VR) (22) 


with /(m) after Table I in (17). The results of using this formula are given for 
two examples in Figs. 2 and 4 (half-body and }-calibre rounded head). The broken 
lines correspond to the constants as proposed by Thwaites“'”, a=0-45 and b=6. 


As he was interested mainly in (plane) cases with increasing pressure—whereas 
in these examples the region with decreasing pressure is of importance—better 
results (single marked points in Figs. 2 and 4) are reached by a=0-45 and b=5:5*. 


Thus it is shown that this simple formula, (18) or (20), can be used for axi- 
symmetric flow as well as for plane flow with an accuracy sufficient for most 
practical purposes. But it seems advisable to make b=5:5 as long as dU/ds>0 
and afterwards b=6 for dU /ds <0, with a=0-45 in both cases. 


APPENDIX 


For practical use some known formulae are given, (a) for the flow near the stagna- 
tion point and, (b) to estimate the indifference point. 


*Note added in proof: Meantime E. TRUCKENBRODT, Ing.-Arch. 20 (1952), p. 211, proposed 
the same method with a=0-47 and b=6-0. 
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STAGNATION POINTS 
= @PLANE FLOW 
: @ AXIALLY SYMMETRICAL FLOW 
e 
\ FLAT PLATE (PLANE FLOW) 
q 
2-84 
\ 


2-04 
J 


Fig. 5. 


Diagram for estimating the neutral stability point from the profile parameter H and the local 
Reynolds number of the boundary layer (for Hartree profiles). 


(a) Flow near the stagnation point 


(a) Two-dimensional : 
U/U,=u,(s/R’), H=2-21,, X=0-0845s, 


0-292 
pU? v Ju, s/R’” R’ VR Vu, 


(8) Axi-symmetric: 


U/U,=u, (s/R’), H=2:29,, 2=0-32,, X=0-0613s, 


0-248 
36 


G 
% 
2: % 
% 
\| 
% 
ge 


METHOD FOR CALCULATING 


LAMINAR BOUNDARY 


LAYERS 


(b) Neutral stability point 


Schlichting’ and Pretsch“*) have shown that it is possible to attach to each 
boundary layer profile approximately a critical Reynolds number : 


Racrit =(U 


If the actual R, is smaller, the laminar layer is stable in any case; if R,>R,,r, the 
layer is unstable to certain disturbances so that transition can occur. The real transition 
to the turbulent layer will take place somewhere between this neutral stability point and 
the theoretical separation point of the laminar boundary layer. 


Usually R.,,rit is given as a function of the parameter A. Yet as A is determined by 
H for one parametric classes, Fig. 5 gives directly R.,,.;;—actually the logarithm of it— 
depending on H. 


This diagram holds for axially symmetrical flow as well as for plane flow. It is to 
be seen that a good approximation is given by the straight line corresponding to 


When X (s) and H (s) have been calculated the critical Reynolds number of the main 
flow Rerit=(U, R’/v)cri, Cam be estimated by (24) for any distance s of the neutral 
stability point from the stagnation point, respectively the position of the neutral stability 
point s for a given Reynolds number R. For this it is practical to write (24) as follows 


log,» =7°55 6-95 (H 2:2)- log, i « 
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Buckling of Oblique Plates with Clamped 
Edges Under Uniform Shear 


W. H. WITTRICK 


(Department of Aeronautical Engineering, University of Sydney) 


Summary: In a previous paper") the author derived, by the Rayleigh-Ritz method, 
a determinantal equation for obtaining the critical magnitude of any uniform system 
of edge stress applied to a clamped oblique plate. In the present paper this equation 
is used to derive values for the buckling stresses of clamped oblique plates subjected 
to pure shear along and perpendicular to two edges of the plate. 


As would be expected, it is found that reversal of the direction of the shear produces 
a change in the critical value. The lower value occurs when the shear is tending to 
increase the obliquity of the plate and critical values corresponding to this direction are 
calculated for 45° oblique plates with side ratios varying between 3/5 and 5/3. The 
critical values obtained for the reverse direction were obviously very inaccurate, due to 
the inadequacy of the series which was used to represent the buckling mode. As yet 
no satisfactory method has been found for overcoming this difficulty. 


It is shown that for 45° oblique plates with clamped edges the maximum ratio 
between the larger and smaller critical shear stresses is 2-69, as compared with 8-00 
when the edges are simply supported. These are the limiting values for very short plates. 


1. Introduction 

In a previous paper”) the author considered the stability of a flat oblique 
plate (i.e. in the shape of a parallelogram) with all edges clamped, subjected to any 
uniform distribution of edge traction. The Rayleigh-Ritz method was used and 
a determinantal equation, giving the critical combinations of edge traction, was 
obtained. This was used to obtain numerical results for the critical stress when 
the plate was subjected to uniform compression in a direction parallel to two of 
the sides. 


In the present paper the analysis of the previous paper is used to obtain some 
numerical results for the critical shear stresses of clamped oblique plates, when the 
shear is applied along and perpendicular to two of the sides of the parallelogram. 
This means that along these two sides the traction is a pure shear stress, whereas 
along the other two inclined sides the traction consists of both shear and direct 
stresses of such magnitude that every infinitesimal rectangular element is in a 
state of pure shear (Fig. 1). The computation is confined to plates having an angle 
of obliquity of 45° and side ratios varying between 3/5 and 5/3. 


Unlike a rectangular plate, it was found that a reversal of the direction of the 
Shear stress causes a change in its critical value. The lower value occurs when the 
Shear is tending to increase the obliquity of the plate. Numerical values were 
obtained for the critical shear stresses when applied in this direction, and involved 
the calculation of the roots of eighth-order determinants. 


Received June 1953. 
[The Aeronautical Quarterly, Vol. V, May 1954] 


ent 
ary 

of 

on 

d.) 

re 

ar 

al- 

ry 

| 

39 


W. H. WITTRICK 


Details of oblique plate. 


The critical values obtained for the reverse direction were obviously very 
inaccurate; this could be seen by an inspection of the convergence of the roots as 
the order of the determinant was increased. In an attempt to discover what order 
of determinant was necessary to obtain an accurate value, the case of an equal-sided 
oblique plate was treated more fully. The determinant was extended to twelfth 
order, but the roots still showed little sign of converging. It was therefore decided 
that the series of Iguchi functions which was used to express the mode was 
inadequate unless a very large number of terms was used. Accordingly, an 
alternative series was tried, but this gave even worse results. 


At the moment therefore no satisfactory results for the larger critical shear 
stresses have been obtained. It is clear, however, that as the length of the plate is 
increased indefinitely, so that an infinite strip under edge shear is obtained, the two 
critical shear stresses tend asymptotically to the same value. It is also to be 
expected that the ratio between the two shear stresses will increase as the plate 
becomes shorter, the limiting case being an infinitely wide oblique plate. For an 
obliquity of 45° this case can be readily solved and the results show that, when 
all edges are clamped, the limiting ratio between the two critical shear stresses is 
2°69. The corresponding figure for a plate with all edges simply supported is shown 
to be 8. It therefore appears that the phenomenon is likely to be more marked in 
a simply supported plate than a clamped one. 


It is shown in Section 3 that on the whole the smaller critical shear stress of 
an oblique panel is somewhat less than that of an equivalent rectangular panel 
of the same weight. A comparison between oblique and rectangular panels in 
compression shows exactly the reverse effect. 


Notation 
a,b lengths of the sides of the parallelogram 
a angle of obliquity of the parallelogram 
t thickness of the plate 
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flexural rigidity of the plate 
=al/b 

Cartesian co-ordinates (Fig. 1) 
oblique co-ordinates (Fig. 1) 


non-dimensional oblique co-ordinates equal to u/a and w/b 
respectively 


Cartesian co-ordinates (Fig. 6) 
=X /(acos Y/(acos 2) respectively 
lateral deflection of the buckled plate 
critical shear stress (Fig. 1) 
=f, cos* «/D 
(acos 
Ty 
ratio between the larger and smaller critical shear stresses of an 
oblique plate 


a,, coefficient of the general term in the series representation of 
w (equation (1)) 


F,,G, functions defined by equations (3) or (4) 
Ampnq the general element of the stability determinant 
Imp” definite integrals defined by equations (7) 


m,n,p,q positive integers 


The Stability Determinant 


The plate under consideration is shown in Fig. 1; the angle of obliquity is « 
and the lengths of the sides are a and b. It is subjected to a pure shear stress, f,y, 
which is considered positive when in the direction shown. Oblique co-ordinates 
u and v are used, and the non-dimensional oblique co-ordinates, € and », are 
defined by the equations 

£=u/a, n=v/b. 


The equations of the sides of the plate are then £=0, 1 and »=0,1. The boundary 
conditions are 

w=0w/0&=0 

w=0w/dn=0 


where w (é,7) is the lateral deflection of the plate when buckling occurs. 


A suitable expression for w is 


w= 
p= 


OG.) 
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where, in order to satisfy the boundary conditions term by term, the set of functions 
F and G are such that 


F, (0) =F, (1) =F,’ 0) =F,’ (1) =0 


G, (0)=G, (1)=G,’ (0)=G,' (1)=0 


A possible form for the F and G functions is that first used by Iguchi’, 


IP +(— IP pre | 
Gu =n sin gry | 


These were used by the author") for determining the critical compressive stresses 
of oblique plates. 


Alternatively we may take 


F,(€)=sin xésin (4) 
G, (y)=sin sin gry f ve 


These functions lead to a somewhat simpler stability determinant by virtue of their 
near orthogonality. It will be seen later, however, that for a given order of 
determinant they give less accurate critical stresses. This was also found to be 
true for the case of an oblique plate under uniform compression. 


The condition for neutral stability of the buckled plate now leads to the 
determinantal equation"). 


Ampna | 


m,n 


where 
+A? Ting” +2 (1 +2 sin? &) Imp? Ing? — 4 sin - 


+ r + 2 Bun (A — sin x (6) 
1 1 
and | = (é) dele dé= Gn (n) dn’ G, () d ] . (7 
0 0 


The constant A is the side ratio a/b, while 8,, is the buckling stress coefficient 
defined by the equation 


2 3 
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Fig. 2. 
Curve showing the smaller critical shear stress of 45° clamped oblique plates. 


The buckling mode was shown" to be either symmetrical or anti-symmetrical 
about the centre of the plate. For a symmetrical buckle we include in the deter- 
minant only those elements in which both (m+n) and (p+q) are even, while for an 
anti-symmetrical buckle we take (m+n) and (p+q) to be odd. The governing mode 
is that which gives the lowest critical stress. 


For given values of A and 2; every element of the determinant in equation (5) 
contains a constant term and one proportional to the coefficient 8,,. The numerically 
smallest positive and negative values of 8,, that make the determinant zero give 
the required critical stresses. In all cases it appeared that the negative critical stress 
is numerically smaller than the positive one. 


3. The Smaller Critical Stress 


The numerically lowest negative roots of the stability determinants were 
calculated by a process of trial and error, as described in Ref. 1. The results for 
45° oblique plates, as obtained from eighth-order determinants, are given in 
Table I. The rows and columns included were those corresponding to m, n, p 
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TABLE I 
CRITICAL STRESSES FOR CLAMPED OBLIQUE PLATES. 2=45°. 


A=a/b 0-6 1:0 135 1-667 


B,,| (m+n) and (p +q)even 11-37 8-99 8-60 8-46 7-61 
™* | Anti-symmetrical buckle. 

(m+n) and (p+ q) odd |  —-7°51 7:16 


and qg<4, and the Iguchi functions (equations (3)) were used to represent the 
buckling mode. 

It will be seen from Fig. 1 that a negative shear stress is one tending to increase 
the obliquity of the plate. The results of Table I are shown in Fig. 2, which gives 
the curve of (G,,/%*)sec plotted against the side ratio A. This form for the 
buckling stress coefficient is chosen for the following reason. It will be seen from 
equation (8) that 

Bua 


(b cos 2) t 


sec =D 


ASYMPTOTE AS 
A> cd 


L 4 i i i. 


L 
2.0 
Seca = a/(b Cose) 
Fig. 3. 
Comparison between the smaller critical shear stresses of 45° clamped oblique 
plates with those of rectangular plates of equal weight. 
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07779 


Fig. 4. 
Infinite strips pnder shear. 


Now bcos z is the width of the plate in a direction perpendicular to the sides of 
length a. Hence for a very long plate (A —> 00) the value of (G,,,/7*)sec z becomes 
independent of « and equal to the value 8-98 for an infinite strip with clamped edges, 
of width bcos, under uniform edge shear’. The curve for 2=45° is therefore 
asymptotic (for large values of A) to the curve) for z=0 which is also shown 
in Fig. 2. 


The volume of an oblique plate is equal to tabcosz. Thus, in order to make 
a direct comparison between the critical shear stresses of oblique and rectangular 
panels of equal weight the buckling stress coefficient (8,,/”) sec x should be plotted 
against the ratio a/(bcos2)=Aseca. This is shown in Fig. 3 and it will be seen 
that on the whole the 45° oblique plate has a somewhat smaller buckling stress than 
the equivalent rectangular plate. 


This result might have beén anticipated in view of the fact that the buckling 
mode of an infinite strip in shear contains nodal lines similar to those shown 
sketched in Fig. 4(a). If the strip is subdivided into oblique panels by a series of 
equidistant stiffeners in the same general direction as the natural nodes of the 
infinite strip (Fig. 4(b)) it might be expected that the critical shear stress would be 
less than if the subdivision were into rectangular panels (Fig. 4(c)). In a similar 
way it might be expected that for a strip subdivided into oblique panels as in 
Fig. 4(d) the critical shear stress would be higher than for the rectangular panels, 
although the results obtained for this case are not sufficiently accurate to confirm 
this conclusively (Section 4). 


It is worthy of note that if an infinite strip is subjected to compression parallel 
to the edges the natural nodal lines are straight and perpendicular to the edges. 
Thus we might expect that in compression an oblique plate would have a greater 
buckling stress than the equivalent rectangular plate. This can be seen from the 
results of Ref. 1 for clamped oblique plates, which are shown plotted in Fig. 5. 
The work of Anderson’ on the compressive buckling of continuous flat sheet 
sub-divided into oblique panels gives further justification for this conclusion. 


It is therefore concluded that, although obliquity gives an increase in the critical 
compressive stress of a panel of given weight, it also causes a decrease on the whole 
of the smaller critica] shear stress. 
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Comparison between the critical compressive stresses of clamped oblique 
and rectangular plates of equal weight. 


4. The Larger Critical Stress 

The lowest positive root of the stability determinant for the symmetrical 
buckling of an equal-sided panel of 45° obliquity (A=1, «=45°) was calculated in 
exactly the same way as the negative root. The Iguchi functions, given by 
equations (3), were again used to express the mode, but this time the roots of 
successive determinants up to the eighth order showed no sign of convergence, 
i.e. each unit increase in the order of the determinant gave a substantial decrease 
in the root. In an attempt to discover what order is required to give an accurate 
value, the determinant was extended to the twelfth order* and the root again 
calculated. The results for the determinants of ninth to twelfth orders are given 
in Table II. 


It will be seen that there is still no sign of convergence of the successive roots, 
and it was therefore concluded that a very large number of terms in the series for 
w is required to represent the mode adequately. 


*This corresponded to the addition of terms in the series, equation (1), involving the coefficients 
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Fig. 6. 
Limit of a 45° oblique plate under shear as A->0. 


An alternative form, given by equations (4), was then used for the F and G 
functions and the roots again calculated. It was found, however, that the accuracy 
was worse than that obtained by the use of the Iguchi functions. 


Thus, at present, no satisfactory results have been obtained for the larger 
critical stress. 


5. The Limiting Values of the Ratio between the Larger 
and Smaller Critical Stresses 


It is obvious that, as A—> 00, so that an infinite strip under edge shear is 
obtained, the direction of shear becomes immaterial. Hence if we define p as the 
ratio between the numerical values of the larger and smaller critical stresses we have 


as 
This result is true whatever the boundary conditions of the panel. 


It is also to be expected that the ratio p will increase as the side ratio A 


‘ decreases, and that p will have its greatest value as A —> 0; we then arrive in effect 


at an infinite strip of width acosz. Considering the case of an oblique plate in 
which «= 45°, it is evident that if f,, is positive this strip is subjected to a longitud- 
inal compressive stress f,, and a lateral tensile stress also equal to f,,, as shown in 


TABLE II 
Order of determinant 9 10 11 12 
Bay 17:70 17-41 16-29 13-03 
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Fig. 6(a). If, however, f,, is negative the strip is subjected to a longitudinal tensile 
stress f,, and a lateral compressive stress f,,, as shown in Fig. 6(b). 


Consider first the case where f,, is positive, as represented by Fig. 6(a). Taking 
axes X and Y, where X is along the centre line of the strip and writing 


X 


Y’= 
acos acos 


the differential equation for the lateral deflection w when buckling occurs is 


son) w= ow 


where k is defined by the equation 


On using equations (8) and (10) we see that 


The buckle will be sinusoidal in the X direction and accordingly we write 
where {2 is a constant defining the wavelength of the buckle in the X direction. 


Substitution of equation (12) into equation (9) yields an ordinary differential 
equation in W whose general solution is 


W=A cosh2 6Y’+Bsinh2 6Y’+Ccos2¢Y’+Dsin2 Y’ 


=F 80% | 
(14) 


G +8 (k+2 a) | 


Suppose now that the edges £=0 and 1 of the oblique plate are clamped’, i.e. 
both edges of the infinite strip are clamped. The boundary conditions for W 
are then 


W=dW/dy’=0 at Y=+}.. . (15) 


*Note that no mention is made of the conditions along the other two edges of the oblique plate. 
These are immaterial when A>0. 
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Fig. 7. 
Curve giving the smaller critical shear stress of short 45° 
clamped cblique plates. 


Also the lowest buckling stress will be associated with a mode symmetrical 
about Y=0 so that B=D=0. Equations (13) and (15) then lead to the 
transcendental equation 


It will be seen from equations (14) that equation (16) defines k as a function of 2. 
It is now a matter of choosing {2 so that k has its minimum value. By trial and 


error it was found that 


occurring at a wavelength defined by 2=2-05 approximately. 


If on the other hand the edges £=0 and 1 of the oblique plate are simply 
supported, the boundary conditions for W become 


W=@W/dY?=0 at Y’= +}. ‘ (18) 


Again the governing mode will be symmetrical, so that B= D=0, and equations 
(13) and (18) now give the equation 
cos ¢=0. 
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It follows therefore that ¢=7/2, and on using the second of equations (14) we 
find that 

(Q? + 1)? 


On minimising k with respect to 2 we obtain 
k=8 
occurring at a wavelength defined by 0= /3. 


Consider now the case where f,, is negative, as represented by Fig. 6(b). It is 
obvious that the lowest buckling stress corresponds to a mode which is cylindrical 
(i.e. w is independent of X), in which case the longitudinal tensile stress has no 
effect. Thus the value of the lateral compressive stress f,, required to cause 
buckling is given by : 


if the edges of the strip are clamped, and by 


if the edges are simply supported. 


The results of this analysis are summarised in Table III. This shows the 
limiting values of (8,,/%?)sec z, as given by equations (11), (17), (19), (20) and (21), 
corresponding to the larger (positive) and smaller (negative) critical stresses of a 
45° oblique plate, and also the ratio p between the two stresses, as the side ratio A 
tends to zero. 


The two values of p given in Table III indicate that for an oblique plate of 
given dimensions there is a considerably greater difference between the two critical 
shear stresses if the edges are simply supported than if they are clamped. 


The results of this analysis also enable the results for the smaller critical stress 
to be extended to values of A less than those covered by Table I. This may be 


TABLE 
CRITICAL SHEAR STRESSES FOR A 45° OBLIQUE PLATE AS A—>0 


oe Edges €=0 and 1 Edges €=0 and 1 
Boundary conditions clamped simply supported 
Langer 10-76 8 
critical 
B stress 
—> sec a 
ad Smaller 4 1 
critical - 
stress x 
p 2-69 8-00 
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done by plotting A? (G.,/%?)sec 2(=k) against A, as shown in Fig. 7. The value 


for 4=0 was taken from Table III, while the values for A=0-6, 0°8 and 1:0 were 
calculated from Table I. The values for A lying between 0 and 0-6 may easily be 
interpolated from the curve. 
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Theory of an Accelerated Slender Supersonic 


Aerofoil 


GEOFFREY L. SEWELL, M.A., Ph.D.* 


(University College, Hull) 


Summary: The problem of a uniformly accelerated slender supersonic aerofoil is 
solved in this paper by a perturbation method for the cases when the acceleration is 


small. 
shock line. 


Notation 
Ox, Oy 


k 


t 
Po> Po» a, 


u(=(u,v) ), p. p 


V 
¢(=V) 
£=€ (0) 

C 

F (x) 
L 

L, 


Full account is taken of the effect of the acceleration on the form of the 


Tig. 1. 


axes of reference, O being the nose of the aerofoil and Ox 
its axis 


unit vector perpendicular to Ox and Oy, making a right- 
handed set with the directions of Ox, Oy, respectively 


time variable 


constant pressure, density and speed of sound, respectively, 
upstream of the shock 


perturbation velocity of flow, pressure and density, 
respectively, downstream of the shock 


velocity of the aerofoil, its direction being opposite to Ox 
uniform acceleration of the aerofoil 

distance travelled by the aerofoil at time ¢ 

given by a,=V sin u 

upper side of the aerofoil profile 

given by y=F (x) as the equation for the form of C 
shock curve 

the line y=x tan p. 


*At present at Marischal College, Aberdeen. 
Paper received March 1953. 
[The Aeronautical Quarterly, Vol, V, May 1954] 
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ACCELERATED SUPERSONIC AEROFOIL 


Theory 


The boundary conditions at the shock’ are given by 
a,?:p, on L; 
and the boundary condition at the aerofoil is 
v=VF’(x) on y=0. 


The linearised equations of flow downstream of the shock are 


div u=0, 
ou ou 
grad p+ p, a + v=) =0, 


0 0 
and (< 
Hence, by equations (1) and (3), 


dp Op 
ax +p,a,* divu=0. 


Moreover, by equation (2), p and u may be expressed in the forms 


oH 0H 


and u=grad H+-curl[ |. 


respectively, and therefore, by equation (4), 


0 
2V72 


Now we know the solution of this equation, subject to the boundary conditions 
(A) and (B), when V is constant"), and this suggests letting V replace ¢ as an 
independent variable, so that 0/dt=<0@/0V; whereupon, neglecting <*, equation (7) 
takes the form 


oH 
2 2 2 


(8) 
If, further, we let H=H,+¢H,., where H, is the solution when V is constant, then 


H,= — VF (x— ycot p) tan p, ‘ (9) 
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and H, is given by the equation 


OH, oH, 
dy? 


Thus, from equations (9) and (10), it may be seen that 


OH, 
— COs? 


sin? ay? = [2 tan’ (x—ycotp)+ 


+2 y sec? uF” (x—ycoty)], 
the general solution of which is 


H, = —3V-' sec? » tan’ (x+ ycot F (x—y cot »)— 
-4V~—'sect ncot . (12) 


and hence 
H=-—VF (x—ycotp) tan (x—-ycot »)+B(x+ ycot sec? » n(x + 
+ ycot F (x—ycot sect cot » (x—ycot p)]. (13) 
The function z, introduced in equation (6), must be of the order of ¢, since 


it is zero when V is constant; and since none of the other functions concerned 
involve t, save in V and p, it follows that « must be a function of y only, so that 


We have now three arbitrary functions, A, B and «,, which we now determine 
from the boundary conditions (A) and (B), the result being that 


a,’ (y)=V—' sec” u F’ (0) y, 
B’ (x+y cot sec? tan® F’ (0) ycot p), 
and A’(x—ycotu)=}V—' sec? » tan*® u[3 (x— y cot ») F’ (x—ycot | 
— 3F (x—ycot (0) (x—ycot »)]. j 
Thus, H and «a are given by the expressions 
H= —VF (x—ycot ») sec? [(x? + cot? F’ (0)- 
(r—y cot p) 
— 6y cot uF (x— y cot F (t) dt —2y? sec? cot* F’ (x— ycot »)], 
and sec? u F’ (0). 
This, then, completely determines the solution, from which it is seen that waves 
are in general propagated along both sets of characteristics, owing to the accelera- 
tion, and that vorticity is produced. 


REFERENCE 


SEWELL, G. L. (1950). Theory of an Oscillating Supersonic Aerofoil. Aeronautical Quarterly, 
May 1950. 


54 


| 


The Delta Wing in a Non-Uniform 
Supersonic Stream* 


G. N. LANCE, M.Sc., Ph.D. 


(Department of Engineering, University of California, Los Angeles) 


Summary: A generalised conical field theory is developed and is applied to delta 
wings in a non-uniform stream. It is shown that a non-uniform stream may be 
characterised by the downwash at all points in space. The lift of a delta wing is found 
when the downwash in the wing plane is given as a power series in the co-ordinates in 
the wing plane. The basis of the conical field theory is described in some detail but the 
results only of the calculation of the lift distribution for various downwashes are given. 
The solutions of certain integral equations, required in the calculations, are given in 
the Appendix. 


1. Introduction 


The co-ordinate system is shown in Fig. 1. The apex of the wing is taken as 
origin, the x-axis is drawn in the direction of the free stream (the wing is placed 
symmetrically with respect to this direction), the y-axis is drawn in the starboard 
direction and the z-axis completes a right-handed co-ordinate system. 


It is shown in Section 2 how non-uniform flows arise in practice and it is also 
shown that the consideration of such flows reduces to the following problem:— 
Given the downwash w over the wing area, assumed in the z=0 plane, as a power 
series in x and y, determine the x-component of velocity; from this the lift of the 
wing can be obtained. Linearised supersonic flow theory is used throughout. 


There are two types of wing to consider: (a) delta wings with subsonic leading 
edges or, in other words, those wings which are swept inside the Mach cone from 
the apex; (b) delta wings with supersonic leading edges, that is wings for which the 
leading edges lie outside the Mach cone from the apex. 


Type (a) are the only wings considered here, because the method of sources and 
sinks can be applied to determine the lift of wings with supersonic leading edges. 
To solve the problem for wings with subsonic leading edges, conical flow theory is 
used. The theory of higher order conical flows':’*':'*) is used, since the well- 
known first order theory is not adequate to deal with the present problem. (An 
exact definition of higher order conical flows is deferred until Section 3.) In this 
paper the general procedure only is described and the final results are listed in 
Section 10. The reader is referred to Ref. 6 for further details of the working in 
any particular case. 


*This paper is a shortened version of Part III of a Ph.D. thesis approved by the University 
of London; it was written while the author was at King’s College, London. 

Paper originally received July 1952. 

[The Aeronautical Quarterly, Vol. V, May 1954] 
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Notation 


t 
k 


I(x, y) 
n 


u’,u”’,v’, etc. 


u, (n) 


G. N. LANCE 


Fig. 1. 
The delta wing and associated notation. 
Mach lines are denoted by — — — — — 
2y is the angle between the leading edges of the wing 
coty=k, and cotu=B. 


parameter, less than 1, related to the angle of the delta wing by 
a=B/k 


constants used in the results of the calculation of the lift distribution 
related to y by coty=k 

non-dimensional lift distribution of the delta wing 

“order” of the conical flow 


components of the perturbation velocity given by u=0¢/0x, 
v=09/dy, w=09/0z 


defined by u’=09'/0x, u” =09” /0x, v’=09'/dy and so on 
value of u(»,¢) on real axis (see equation (24)) 


constant which is the coefficient of x’y* in the expression for w’ 
(see equation (3)) 


rectangular Cartesian co-ordinates (Fig. 1) 

defined by B?=M?-1 

constants given in Section 10 

complete elliptic integral of the second kind, modulus (1-—a’}! 
complete elliptic integral of the first kind, modulus (1-— a’)! 
Mach number 

mean free stream velocity (Section 2) 


analytic functions of complex variable 9 or c, real parts of which 
are u(y, v(y, 0) and 


denotes the quantity (x*— k’y’} 
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y half the angle between the leading edges of the wing 
7, real and imaginary parts of o 
Mach angle 
£,9 co-ordinates defined by (6) 
complex variable defined by if 
@ perturbation velocity potential 
¢’,.¢”,9¢”” velocity potentials defined in Section 2 and satisfying equation (2) 
© complex variable related to « by (21) 
® complex conjugate of 9 
®,4 analytic function of o, with p+q+r=n 
symbol denoting “ principal value integral” 
HR “the real part of” 


2. Non-uniform Flows 


A non-uniform flow is one in which the free stream velocity is a function of 
the space co-ordinates. It is assumed that the perturbation velocity potential and 
velocity components are independent of time. 


Suppose now that we have an aeroplane flying with velocity U into still air, 
then it is simpler to consider the aeroplane stationary with a wind whose speed is 
U blowing past it. The free stream flow past the mainplane can be considered as 
uniform, if the presence of propellers and so on is neglected; the lift on aerofoils 
in such a flow is well known. However, this mainplane will certainly disturb the 
free stream and on reaching the tailplane the wind will be non-uniform. The case 
of a rectangular tailplane, with large aspect ratio, has been solved by Mirels”° and 
we are here considering a delta shaped tailplane. The non-uniformity is best 
described by the downwash, which is produced by the mainplane, at the surface of 
the tailplane and which would be present if the tailplane were absent. The 
methods of calculating the downwash behind a supersonic wing have been given by 
Robinson“, Robinson and Hunter-Tod"*’ and Ward"*’. A non-uniform flow also 
occurs in a wind tunnel which is imperfect, in that the liners do not succeed in 
producing a uniform flow. In a perfect wind tunnel there will be no downwash, so 
again the non-uniformities in the stream may be described by stating the downwash 
at all points in the field. The results obtained in the present paper can be used to 
find the lift of swept wings in supersonic flow when aeroelastic effects are taken 
into account. A few further results are required which can be derived by methods 
similar to those used here. 


The fundamental linearised equation for the velocity potential, 9, is 


2 2 2 


ay? Ox?’ 
where B?=M?—1, and M is the Mach number. 
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Since (1) is a linear partial differential equation, the ¢ may be split up into 
various parts and if each part is a solution of the equation so is their sum. We write 


+9" . . . (2) 
and define the three components in the following way :— 


(i) 9 is the perturbation velocity of the flow, due to the upstream wing, 
when the tailplane is assumed to be absent. 


(ii) ” arises out of a cancellation of w’=09’/dz (the perturbation velocity 
of the non-uniform stream in the z-direction) at each point of the tailplane 
plan form area. In other words ¢” is to be constructed such that 
09” /0z= —09’/0z for all points in the z=0 plane which lie within the 
tailplane plan form area. 


(iii) 9”’ is the velocity potential due to the tailplane, at given angle of 
incidence, in a uniform stream whose velocity is U. 


When we consider the downwash at all points of the tailplane surface due to 
¢, 9%” and ¢” it is clear that the boundary condition—the flow is tangential to the 
tailplane—is satisfied. The values of ¢’ and ¢”’ will be assumed known as these 
have been given elsewhere: '*:**:'*. We shall assume that w’ is known and 
expanded in a series of the form 


W = Woot t+ Wor Wo X? 


where the w,, are constant coefficients. 


The problem thus reduces to finding the value of u” =09”/0x at all points of 
the tailplane, given the value of w”, where w” = — w’ and w’ is given by (3). Having 
found u”, this will immediately give the additional lift produced by the non- 
uniformities. It might be thought that it is sufficient to treat the flow as symmetrical 
about the Oxz plane but this assumption will not be made because, although it is 
justifiable in level flight, when the aeroplane is banked the ailerons on the mainplane 
are brought into action and the flow will lose its symmetry. Furthermore, in the 
wind tunnel and aeroelastic applications of this theory there are no grounds for 
assuming symmetry. 


3. General Discussion of Conical Flows 


The remaining sections are concerned only with the tailplane and with the 
velocity potential 9”, defined in Section 2, so it will be convenient to drop the 
dashes and let 4 be the velocity potential whose downwash w is given in the form 


W = — (Woo + X+ Wo X? XV + +....). (4) 


Furthermore we shall call the aerofoil a wing, instead of a tailplane, and (4) gives 
the downwash on the wing. We also confine our attention to delta wings, 
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symmetrically placed with respect to the free stream direction, whose plan forms 
are isosceles triangles lying inside the Mach cone from the apex. 


First order conical flows have been discussed in very great detail by a number 
of authors, in particular Multhopp"’’, Germain":*), Goldstein and Ward“, 
Stewart"”), Lagerstrom®), Ribner"*’, and Lomax and Heaslet'’. There is little 
point in repeating the majority of the first order theory but we will give sufficient 
to enable the reader to understand the higher order theory without reference to 
other papers. 


We define a conical flow of the n“ order as a flow for which all the n™ partial 
derivatives of the velocity potential ¢ with respect to x, y or z are constant along 
straight lines through a fixed point. With this definition it is seen that the type 
of conical flow, considered by the authors cited, is of the first order, i.e. one in 
which the components of the perturbation velocity are all constant along straight 
lines through a fixed point, the vertex of the delta wing. 


Higher order conical fields have been mentioned by Multhopp"”’, and Lomax 
and Heaslet'’, and Germain”) has also considered such flows but assumed that 
conditions were symmetrical about the plane y=0. This will not be sufficient for 
the present purpose and the theory is developed without any such assumptions. 
Multhopp only discussed first and second order theory although he did suggest how 
the work would proceed after that; his suggestion is however rather misleading. 


4. Transformation of the Velocity Potential Equation— 
First Order Theory 


The velocity components of the flow are given by u=09/0x, v=0@/dy, and 
w=09/0z, which all satisfy the Prandtl-Glauert equation. The equation for u is 
taken as typical. 


The co-ordinates (x, y,z) are changed to (R, £,6) by means of the equations 
x=BReoth€, y=Rcosechécos@é, z=Rcosech €£sin 6. (6) 


If a simple cone field is assumed, u is independent of R and it can be verified that 
(5) is transformed to 


One further transformation is now applied, namely, 


_ cos € sinh € 
Cosh? sin? 6 = cosh? € — sin? 
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This change of variables does not change the form of (7), which becomes 


It can be shown that the transformations have reduced the problem to a two- 
dimensional one in the o-plane. The «=7+i€ plane is such that the Mach cone 
has transformed into those parts of the real axis for which —cO<=7»=<-—1 and 
1<y»<0oo and the delta wing appears as that part of the real axis for which 
|» |<a, where a=B/k. (ais less than 1 for a wing with subsonic leading edges.) 


. B By? B27?) 


Since we are only interested here in the value of u and w on the delta wing, 
which is given by | 7 |=| By/x|<a, (=0, we are only concerned in the (n, ¢) plane 
with the values on the real axis. This fact makes for considerable simplification, 
since we need only evaluate the derivatives 0/0x, 0/dy and 0/0z on the (=0 axis. 


The results, which have been obtained before by Multhopp"”, are 


0 Bo 
a 

and (14) 


all of which apply only when (=0 and {| 7|<1. Remembering that the motion is 
irrotational, the following three equations connect u, v and w:— 


Ow dv ou dw ov _ Ou 


Thus, using (12), (13) and (14) the irrotationality conditions become 
B(1—n?} (17) 


respectively, all of which are valid when y and »>0. 
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5. The “Compatibility Conditions ” 


Since, by definition, ¢ is a homogeneous function of x, y and z of degree one, 


in First Order Theory 


x2 =XU+ yu+Zzw, 


Ox oy 0z 
and therefore pre + ydv + zdw=0, 
or Bcosh du+ cos dv +sin 6 dw=0. (19) 


From the fact that the velocity components satisfy Laplace’s equation (7) inside the 
Mach cone it follows that u (0, €), v (6, £) and w(@,€) are the real parts of analytic 
functions U (9), V (Q) and W (0), where O=6+i€. 


Equation (19) can now be written 


Bcosh € dU +cos 6 dV + sin 6 dW =0, 


and if © denotes 6—i€, the complex conjugate of 9, it follows that 


os5 | Boos dU +cos dV +sin aw | + 


+sin BsinS au sin dV +cos dW =(Q. 


Each bracket must vanish independently, and so 


dV dw 
. . . (20) 


These are the compatibility conditions and they connect the expressions for u, v 
and w. It will be convenient to apply the further transformation (8), namely 


In this case (20) becomes 
ir dW 
—-BdU=cdV=- 


When considering the higher order conical fields, more and more analytic functions 
will be introduced and also more and more relationships connecting them. For 
instance in the present case, i.e. first order case, there are three functions, in the 
second there are six, in the third order there are ten and in the n™ order there are 
4(n+1)(n+2) analytic functions. For this reason we introduce a notation which 
can be extended to n™ order flows. The equation (22) will be written 

io d® 


‘ 
)) 
_| 
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where the ®’s are analytic functions of « and where we have replaced dU by d®,,,, 
dV by d®,,, and dW by d®,,,. This notation can be extended to the general case 
(see equation (50)). 


It is satisfying to note that the three relationships (16), (17) and (18) are in fact 
included in the compatibility conditions (22), bearing in mind that o=n+i¢ and 
that u= & U and so on. The remaining equations of (22) merely give connections 
between the imaginary parts of U, V and W which are of no interest to us. 


6. Expressions for Certain Derivatives on the ‘=0 Axis 


It is as well to pause here and summarise the problem. We wish to find the 
lift distribution on a symmetrical delta wing, which is swept inside the Mach cone, 
being given the downwash over the wing as a power series in x and y of the form (4). 
We are, therefore, interested in particular in the relations (22) which connect the 
downwash w with the lift distribution (which is a linear function of u). Examining 
(22) it is seen that 0/0¢ appears in this equation, and since w is given in the 
z=0(¢=0) plane we shall require an expression for 0 (u,v, w)/0¢ on (=0; the next 
task is to obtain an expression for 0 (u,v, w)/0¢ on (=0 in terms of the value of 
(u,v,w) on (=0. Again u is taken as typical of the three velocity components. 


For lifting problems 
u(y, +0)= —u(y, —O)=u,(n), (Say), A) 


for |»|<a, where +a are the extremities of the wing which appears as a cut 
along (=0 in the o-plane. Clearly, 


a 


u(y. = (25) 


+0 


gives the boundary values (24) along the wing cut and also it satisfies Laplace’s 
equation (9). It now follows immediately that 


du(n,0) 1 u, dr’ 
a 


—a 


Il 


(26) 


This integral is improper if | |< a and the principal value of the integral is 
defined by 


a 


€ 


lim 


—a 


nte 


Mangler”? has written an excellent report on improper integrals in which he proves 
that all that is necessary is to evaluate the integral in the ordinary straightforward 
way and any infinities which may occur can be rejected. A note on improper 
integrals also appears in the Appendix of Ref. 7. An integral of this sort, in which 


all infinities are to be neglected, is denoted here by the symbol § . 
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Subject to the condition that u,(—a)=u,(+a)=0, equation (26) may be 
written in the alternative form 


ou 1 {du, dy 
dn’ (n—1)’ 


A further relaticn, which can be proved by integrating by parts, and which will be 
used later is 


d u, (n’) dn’ du,(’) dn’ u,(—a)_ u,(a) 


dy 


7. Higher Order Conical Flows 


The theory considered so far will only enable us to solve the problem of 
finding the lift when w= — w,,=constant, on the wing; we now turn to the general 
problem and show how this may be dealt with. By definition, for a conical flow 
of order n the velocity potential ¢ (x, y, z) is a homogeneous function of degree n in 
(x,y,z). It follows that the n™ partial derivatives of ¢ with respect to x,y or Z, or 
any combination of all three, are all homogeneous of order zero; in other words, if 


with p+q+r=n, then all the possible ¢,,, are functions of y/x, z/x and satisfy the 
two-dimensional Laplace equation (9) in the (y,0) plane. As in Section 5 it follows 
that %,,, may be taken as the real part of an analytic function of the complex 
variable o, say, Pyq, (c). 


By Euler’s Theorem on homogeneous functions it follows that 
! = — 


where the bracket is to be expanded algebraically; the equation implies that if all 
the partial derivatives of ¢, namely ¢,,,, are known, the ¢ itself is determined. Some 
particular partial derivatives are denoted by special symbols, for reasons which 
will be apparent later. Let 


for all OS q<n-1, BL 
also let for all OSq<n-l. . 


If u and w are the x- and z-components of the velocity, respectively, then 
n—1 
q=0 
a in the z=0 plane. 
~1)!w= -a-1 
and (n—1)!w (34) 
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These definitions are consistent with the notation in the first order theory, because 
in the first order theory n=1, so gq=0 only; thus u,,=%,5,=0¢/0x, so (33) shows 
that u=u,,=0¢/dx, as it should; also w,,=%.,=09/0z, so (34) gives 
w=W,,=09/0z, again in agreement with the previous work. 


8. The Compatibility Conditions in Higher Order Theory 


The proof of formulae (22) is quite straightforward for first order conical flows, 
but it is not the best proof in extending the methods to second and higher order 
theories. For this reason an alternative derivation of (22) is given, incorporating the 
notation used in Section 7. 


Beginning, as in Section 5, with the equation 
xdu+ydv+zdw=0, . ‘ (35) 


let ®,,,(o)=u+iu,, P,,,=v+iv, and ®,,,=w+iw,, where «=sec®, as before. 
The component of the gradient of u, in any direction is equal to the component of 
the gradient of u in a direction at right angles; similarly for v and w; hence (35) 
implies xdu, + ydv,+zdw,=0, and so 


the dash denoting differentiation with respect to o. Differentiation of (36) with 
respect to x, y and z in turn leads to 


and two similar equations. Thus 


ax 


, , ‘ , 
100 ® 010 ® 


(38) 
To find expressions for d0/dx, dc /dy and 00/0z, we pick out the coefficients of 


dx, dy and dz from the expression for do in terms of x, y and z. Since c=sec® 
it can be shown that, from the transformations introduced by equations (6) and (8), 


Thus (38) shows that 


d® as ix d®,,,, 


~ do 


which is the relation previously obtained (see equation (22)). 
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We now obtain the conditions when the flow is a second order conical flow. 
In this case ¢ is homogeneous of the second degree in (x, y,z). Hence 


4=2! 


and y = Lh. (43 ) 
Therefore 


(44) 


with all the second partial derivatives of » with respect to x, y and z of degree zero. 
Each of these is thus the real part of an analytic function of «. Letting 


= 
ox” * ay 


be the real parts of the six functions 
respectively, (44) becomes 


Differentiating this equation twice, with respect to x, y and z in all possible ways 
gives the following six equations :— 


02a 
pt+qt+r=2; 


so that 


(47) 


By picking out the coefficients of d?x and so on, from do, we can obtain 0°0 /0x? 
etc., and find that the resulting equations can be written compactly by saying that 


all the quantities 
io », (7) 
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with 0<= p+q<2 are equal. It is easily seen that this expression can be generalised 
and, for n order conical flows, we have the very important and useful result that 
all the expressions 


with O<p+q<n are equal. 


We now have all the tools at our disposal and can proceed to solve the problem. 
First order theory deals with the case when ¢ is homogeneous of the first degree; 
so w, the downwash, is a constant on the wing. Second order theory deals with the 
case when w, on the wing, is proportional to x or y; third order theory when w, on the 
wing, is proportional to x’, xy and y’ and fourth order theory will give solutions to ; 
the case when w is proportional to x°, x*y, xy’ and y’, on the wing. The results 
of all these are given in Section 10; as the results indicate, the working is formidable 
even to third degree functions of x and y. 


9. Method of Solution in the General Case 


The working by which the results given in Section 10 are obtained will not be 
given for all cases but the third order case will be used to illustrate the procedure. 

The third order theory will solve the cases in which the downwash is of 
the form 


The unknown function in this case can be taken to be 0?u/0x?, which is homo- 
geneous of degree zero; this function is written %,,,. Thus to solve the case 


W= we have 
in on . . Gf 

and on (=0; 
Do, =0). (53) 


and we require to find ®,,, in |7|<aon (=0. From (49) we pick out the three i 
equations which connect ®,,, with ®,,,, ®,,, and ®,,,. They are 


d®,,, io d®,,,, —io? d®,,, d®,,, (54) 


Equating real and imaginary parts of these equations, and equating the imaginary 
parts Of Poo, and ®,,, to zero (the real parts are $500, $201» $111 ANd 
respectively) we find 
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The remaining equation is not required. 
We see from (51) that 

00 : 

in |n|<a on . (56) 


and from (55), since ¢.,, is zero at y=1 and —2w., at n=a, 


1 1 


Using (27) it can be shown that 


u, ) dn 


so by (56) we require the form of u, (n’) such that 


Cu, u, dn’ 


= constant + x constant in || <a. 


A reference to the Appendix, result (ii), shows that this equation has a solution 
which is either symmetric or anti-symmetric. Only the symmetric solution is 
considered here, namely 


o+Con? 


where c, and c, are constants to be determined. From (58) and (59) we obtain 


Substituting this value of 34,,,/0¢ in (55) gives 


1 
c 
[Got 
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These integrals can be evaluated in terms of the complete elliptic integrals of the 
first and second kinds of modulus (1—a’)!. Equations (61) and (62) become 


Cy (a? — 3) K +2 (2—a’) E] —c, [2a°K —(1 +a”) E]=2w,, . 
and c, (2a* —1) K + (4 7a’ + a*) E]+ c, [a* — 3) K+ 2 (3 5a’ +.a*) E]=0, (64) 
respectively. The solutions of (63) and (64) are 

— 2w,,a’ [a? (3 — K —2 (3 5a? + a*) E] 
[Sa*K? — 8a? (1 + a?) KE — (4 — 19a? + 4a*) 


+ 2w,, (1 — 2a’) K —(4—7a* + a*) E} 
[Sa*K? — 8a? (1 + a?) KE— (4— 19a? + 4a*) 


Co= 


(65) 


Since u=x?u,(y) and the non-dimensional lift distribution /(n) is given by 


2 
= 


the final result for the lift distribution is 


4k (c,x* + c,B*y”) x 


I(x, y)= UB? (x? — k?y?}? (66) 


where c, and c, are given by (65). 


In the following section the remaining results are listed, all working being 
omitted. It is realised that the working left out in the present section is considerable 
but it is given, in full, in Ref. 6. 


10. List of Lift Distributions for Various Downwashes 
given on the Wing 
(1) w=—Wo; 
UkE (x? — k*y?)}! UkEX 


W= — 
4w,, (1 —a’) (2x? k*y’) 


—4w,, (1 — a’) xy 
Uk —(2—a’) E] X° 


I(x, y)= 


2 


— 2w,,a* [a? (3 — 5a?) K — 2 (3 — 5a’ + a*) E] 
[Sa*K? — 8a? (1 + a”) KE — (4— 19a? + 4a*) 


[2a? (1 — K (4—7a? + a) E] 
— 8a? (1 + a) KE — (4— 19a’ 4a) 
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(5) 
4ky (c,x?+¢,B’y’) 


y)= UBX 


—2w, (1 —a?)? 

~ Bla ) 


C3 


4kx (c,x?+¢,B’y? 


+ 2w,.a* [2a?K — (1 +42) E] 
~ B? [5a*K? — 8a? (1+ a?) KE — (4— 19a? + 4a") 


— 2w,.a@* [(3 — a*) K—2 (2—a’) E] 


C4 


[5a*K? — 8a? (1+ a?) KE—(4— 19a? 


(7) w=- 
_ 4k +c,,Bty*) 
UB?X 


I(x, y) 


where 
C, = 8a? (1 — a’) w,, [a® (3 — 9a* — K — (6— 15a’? + Sa* — 4a°) E]/D, 
C,\)= —2(1 — a’) w,, [a? (9 — 45a? — 4a*) K — (18 — 51a® + a* — 8a*) E]/D, 
—2(1 —a’) w,, [a? 9 — a’) K— (3 + 7a? — 2a*) E]/D, 
and 
D, =a (27 — 31a? — 8a*) K? — 2a? (24 — Sa? — 15a* — 16a*) KE — 
—(12—119a? + 151a* — 64a* + 32a*) E?. 


(8) w=—w,,x’y; 


4ky (c,x* +.c,B?xy’ 
y= xy) 


where 

C;= —2w,,a’ (1 — a’) [a? (12 — 31a? + 27a*) K — (24 — 65a’ + 55a* — 6a°) E]/ BD, 
+2w,, (1 — a’) [a? (8 — 21a? + 21a") K — (16 — 44a? + 39a* — 3a*) E]/ BD. 
and 

D,=a* (8 + 31a? — 27a*) K — 2a? (16 + 15a? + Sa* — 24a°) KE+ 


+ (32 — 64a? + 151a* — 119a* + 12a*) E? 
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(9) w=—Wy,pxy’; 
Non-dimensional lift distribution given by (76) but in this case 
C, = 8a* (1 — a’) w,, [a? (9 — a”) K— (3 + 7a? — 2a*) E]/3B?D, 
C19 = — 2a’ (1 — a’) w,, (81 — 65a? + 24a*) K + (3 — 131a® + 136a* — 
— 48a°) E]/3B’D, 


—a’)w,, [a? (27 — + 12a*) K+ 
+ (6 — 55a’ + 65a* — 24a*) E]/3B°D, 


and D, has the same value as before. 
(10) w=-w,,y*; 
Lift distribution is given by (78) but in the present case 


Cc; = 6w,,a* (1 — a’) (1 — 9a”) K — (2 — 7a? — 3a*) E]/ B°D, 
c,= — 6w,,a° (1 — a’) [a’ (9 — a’) K— (3 + 7a? — 2a*) E]/B°D, 


(81) 


and D, has the same value as before. In all these results the modulus of the com- 
plete elliptic functions K and E is (1 —a’)}. 


The results (1) to (4), (6) (7) and (9) have been given before by various authors; 
the present results all agree with these authors. However the complexity of the 
results increases rapidly and it is desirable to be able to check them by an entirely 
independent method. This has been done by the present author using hyperboloido- 
conal co-ordinates; the method is too long to describe here but it is given in full in 
Ref. 6. The lift distributions so obtained do, in fact, completely confirm equations 


(67) to (81). 


APPENDIX 


SOLUTION OF CERTAIN INTEGRAL EQUATIONS 
In Section 9 we had to solve integral equations which were of the form 


=f (n)=a+ Bn+ yn? ; (82) 


Magnus shows how to solve such an equation; we write »=-—acos@ and 
n'=-—acosy¥. If F(¢)=f(—acos¢) may be expanded in a series of the form 


sin  F (¢)=—4 by sinng, 
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where b, is an arbitrary constant*. In the particular case when f(n) has the form 
given by (82), 
F (¢)=a— Bacos + ya’ cos? — da* cos* 
=[a sin ¢—4 Basin 29+ 4 ya? (sin + sin 39) — (2 sin 2¢ + sin 49)]/sin 
n=1 sin p 
Hence the constants b, are 
b,=-2(a+} ya"), b,=(Ba+45a*), b,=-4ya’, 
and b,=0 if r>4. -Thus, from (83), 
—n”)-* [4 b, — 22+ 4 ya") cos + (Ba+ 8a*) cos — 
4 ya’ cos +4 cos 
=a (a? — [4 b,’ — cos 
+ (2Ba — 8a?) cos? — cos* + 28a* cos* 
= (a? — (Cot + +4"), 
where c, is an arbitrary constant and c,=(22— ya’), c,=(28 — 8a), c,=2y and c,=28 


Thus we have the following pairs of relationships :— 


(i) When f(y) =e, 
@ 
(ii) When f(y) =2+ 

(iii) When f(y) =2+By+ yn’, 

+ + 059") (a? — 

(iv) When f(y) yn’ 

8 =(Co + + (@ — 97) 


In all cases if the constants «, 8, y and 4 are arbitrary, as they are in the present 
applications, then all the c,’s are also arbitrary. From these expressions, and using 
the symmetry or anti-symmetry of the required solutions, the solution of the 
integral equation used in Section 9 is as stated. 


A further method of solving (82) is given by Sohngen"*’. Sohngen proves that if 


an J (n’—n) =f ( 


2 
then g(n (a’ — ) = (a?—1)! dn. (84) 


*If f(n) is discontinuous at 7=0, then (83) is not the —— solution. However, the f (7)’s 
considered in this paper are continuous at 7=0. 
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It can be verified that when f(y) has the form given in (82) the expression for 
g(n’) given by (84) agrees with that obtained by using (83). When f() has a form 
in which all the b,’s of (83) are non-zero then (84) is probably the simpler expression 
to use. Such forms of f(y) have not however arisen in this paper. 
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Some Examples of Three-Dimensional Effects 
in Boundary Layer Flow 


J. WILKINSON, M.A. 
(University of Bristol) 


Summary: The three-dimensional boundary layer flow defined by the external 
irrotational velocity components U =  (€), V=an, where 2 is a small positive parameter, 
is investigated with the aid of the boundary layer equations of Howarth’. When 
x ()=€" a solution exact to the first power in z is found. A Pohlhausen method is 
then developed for any function x (£) and applied to the cases in which y (€)=£" and 


x )=1-§ 
Introduction 


The boundary layer flow past an infinite yawed cylinder in a stream has already 
been extensively studied. Such a flow has external irrotational velocity components 
U= x (€) perpendicular to the generators of the cylinder and V =constant parallel 
to the generators, and is one of the simplest cases of three-dimensional boundary 
layer flow. In order to obtain more information about three-dimensional effects 
in boundary layer flow, and especially on separation, consideration is given in this 
paper to another simple boundary layer flow: that defined by the external irrota- 


tional velocity components Ui= ,(£), V=an, where €,» are suitably chosen 
orthogonal curvilinear co-ordinates in the boundary surface and < is a small positive 
parameter. This flow, although chosen chiefly because its three-dimensional effects 
can be readily investigated, is not without practical interest, for it occurs in the 
flow of a stream past a body with a slightly curved leading edge, near the point at 
which the tangent to the edge is perpendicular to the stream. 


In investigating this flow the three-dimensional boundary layer equations of 
Howarth" have been used and the solution expanded in powers of the parameter 
2, terms in a? and higher powers of z being neglected. If 2 is small enough this 
approximate solution will show adequately the trends of the three-dimensional 
effects on the flow. In general it is not possible to discuss a three-dimensional 
boundary layer flow without a knowledge of the geometry of the surface over which 
the flow takes place, for the curvature terms in Howarth’s boundary layer equations 
depend on the shape of the surface. In this investigation these curvature terms are 
omitted, so that the equations apply strictly only to flow past a plane or a cylinder 
with velocities induced by external means. In practice, with many flows of the 
kind considered, the curvature terms are negligible, at least over a limited region. 


In Part I a solution exact to the first power in z is obtained for various values 
of m when x (€)=€". In Part II a Pohlhausen method is developed for any function 
x (€) and is then applied to the flow considered in Part I to give values for the skin 
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friction in reasonable agreement with those obtained from the exact method. 
Finally the Pohlhausen method is applied to the case where x(£)=1-—€. In this 
case it is known from two-dimensional applications that the Pohlhausen method 
gives poor correspondence with exact methods and the results obtained must be 
considered as qualitative. 


In both cases, x(6)=&" and x(é)=1-—€, the skin friction is found to be 
increased as the result of a cross flow of the type assumed. Because the method 
of solution breaks down for the critical value of m when ,(€)=€" and the 
Pohlhausen solution gives increasingly inaccurate results with increase in € when 
x(=1-6, it is not possible to draw any certain conclusion about the separation 
of the flow: but in view of the increase in skin friction it seems reasonable to 
suppose that separation is delayed by the cross flow. , 


Notation 
orthogonal curvilinear co-ordinates 
velocity components inside the boundary layer 
velocity components outside the boundary layer 
coefficient of viscosity 
kinematic viscosity 
a small positive parameter 
parameter in the expression x (£)=&” 
=2m/(m+1) 
boundary layer thicknesses 
=1/q=8/8 
fluid density 


Ww 
al 
b 
ti 
al 
it 
b 
PART I 
U=€", V=an. 
The equations to be considered are 
ov dv 
ou 4. ov + ow =0 
On 
14 


THREE-DIMENSIONAL EFFECTS IN BOUNDARY LAYER FLOW 


where v—->an as (—->+0 
u=v=0 when (=0 


and 1, v/an—1 as €—->+0. 


These equations are derived from the boundary layer equations of Howarth"? 
by taking h,=h,=1 and substituting the pressure term derived from the irrota- 
tional flow U=£",V=an. The assumption that h,=h,=1 means that the equations 
are strictly applicable only to flow over a plane or a cylinder with the given velocities 
induced by external means. 


It is assumed that there is a solution u=u(€,0, v=anv* (€, 0), the assumption 
being justified a posteriori. In this case w=w(é,0). Expanding in powers of a 


u=u, (€, 0) + au, (€ 9 +0 (27), 
v*¥ =v, (€,0) +O(@), 
w=w, 0+ aw, (€,0 +O (2°). 


Substituting and equating coefficients of powers of ~, the equations determining 
Uo, Uy, V, are 


07u, 
av? 


Ou, 


* 


dv, 


=0. 


(y), u,=€F,’(), 


equations (1) and (2) reduce to 


du, Ow, 
where u,=u,=v,=0 when [(=0, 
_ 2m 
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+ FF,” 
equation (3) becomes 
G,’’+F,G,” =0, 
and equations (4) and (5) become 


where F,=F,’=G,=G,'=F ,=F,’=0 when y=0 


and F,— 1, G,—->1, as 


Equation (6) is the Falkner-Skan equation in the form in which Hartree?) 
used it. Hartree has given the solution for many values of 8. Cooke‘) has con- 
sidered equation (7), which can be solved by quadrature, and he gives the solution 
for some values of £. 


The new equation, equation (8), has been solved for a few values of 8 and the 
quantities F,’(y), G,’ (y), F,” (0), G,” (0) are shown in Table I. 


m—1 


so that F,” (0) and G,” (0) are a measure of the change in skin friction caused by 
the three-dimensional motion. 


The values of F,’ (y), G,’ (y), F,” (0), G,” (0), for B=1 have been taken from a 
paper by Howarth’, who calculated them in finding the solution near a stagnation 
point. For 8=—0-18, equation (8) was solved by the Taylor series method and 
for the other values of 8 the central difference method was used. 


When 8=0 the solution is trivial, with 


G,(y)=F,() and 
so that G,’Q)=F,’(y) and F,’(y)=4yF,” (y). 
In this case F,” (0)=0°2348 and G,” (0)=0-4696. 


When 8=2 the solution of equation (8) is identically zero and when 
8=—0-1988, which corresponds to two-dimensional separation, it seems to be 
impossible to find a solution which satisfies the boundary conditions. As 
B—>—0-1988 from above, F,” (0) —> 


As < is positive, and supposed to be sufficiently small for terms in 2? to be 
ignored, the skin friction is increased for each £, the increase becoming greater as 
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TABLE I 
TABLE OF F,’(y) AND G,’(y) 


— 


6905-0 606¢-0 (0) 
9IST-0 8S9L-0 (0) 


(0) ,'9 ANV (0),,'y dO ATAVL 


IN BOUNDARY LAYER FLOW 


0-9 
9-¢ 
(aS 
0-S 
9-4 
ve 
ce 
8-7 
9-7 
0-7 
8-1 
9-1 
0-1 
8-0 
9-0 
v-0 


THREE-DIMENSIONAL EFFECTS 


(4),'9 anv (4),'4 JO 
I 


| 
COD 
CONN AO 
RRARR— 
ow 
~o 
COSCO 
on 
A 
oo 
SSSSSSSS80 
= oO OF 
RRARRAD— 
| 3S 
moO t™ oo 
CO 
| 
| 
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B decreases. The failure of the method of solution for 8= —0-1988 means that 
nothing definite can be said about the effect on separation, but the general increase 
in skin friction, together with the form of failure at B= —0-1988, makes it reason- 
able to suppose that in fact separation will be delayed by a cross flow of this form. 


One of the simplest applications of these results is to the flow near a stagnation 
point on a slightly curved, blunt leading edge—for example, to the flow near the 
forward stagnation point on a thin planetary ellipsoid of revolution placed in a 
stream parallel to an equatorial diameter. Then, with a suitable choice of 
co-ordinates, the external potential flow will be U=£, V=an—that is to say, in 
this case, m=8=1. Using the figures in Table I and those given by Hartree", it 
follows that the skin friction is pv!£{1-2326+0-0478 «+ O (2”)}. 


The fractional increase, due to the introduction of the small cross flow V by 
the curvature of the leading edge, is to the first order in 2, only 0-0388 a. In view 
of the fact that by hypothesis « is small, this change from the value for the two- 
dimensional flow past a straight edge is very small indeed. 


PART II 


With h,=h,=1 in the boundary layer equations, the appropriate momentum 
integral equations are 
3 3 


where 6 and 8 are the boundary layer thicknesses corresponding to the u and v 
velocities respectively. 


Write u= xf (x,A), v=ang(y) 


-(2x— 2x —x) t<8 
2y*+y*) (<v 


where f (x, A) 


x=C/8, y=6/8, A= ‘8? 


If it is assumed that 6 and & are functions of £ but not of y—an assumption which 
can be justified a posteriori by finding such a solution satisfying the momentum 
integral equations—then, expanding in powers of 2, 
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where ¢,, ¢,’, €, are functions of only. Substituting in equations (9) and (10) and 
equating coefficients of powers of « 


1 1 


0 
1 


1} g (px) dx=xy’e,+ 


where p=1/q=8/¥ and = A,=y’.,?/v. 
Equation (11) is the ordinary Pohlhausen equation for the two-dimensional flow 
U= x 

Equation (12) can be written 


d 1 da, 
13 (5, ~ x 


1 


where G(q)=4q f ay. A.) [g — 1] dy 


p+ 
- 


10p 


1 1 1 


if q<l. 
On writing ¢,=k()¢,, equation (13) becomes 


(15) 
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315 * 315 * 9072 


630 2520 4536 


F, (A,)= 


F,Q,)= 


+ iw”) +7 + 3004? ) if 


67 5, 7 4) _ | 
=4(-7 10 364 10 


of 13,201, 7 


The method of solution for any given function x () is as follows: — 


if 


(i) A, is found from the equation 


20,)+% 2h (dy) +2, 


where g(A,), A(A,) are the functions defined in Ref. 5, page 160. 


(ii) p or q is found from equation (14) which, when A, is known, is a first 
order, first degree, non-linear differential equation. As p, q=1 usually, iteration 
is a suitable method of solution, beginning with p(¢)=1. 


(iii) k is found from equation (15) which, after the known values of p and 
A, have been substituted, is a first order, linear differential equation. 


It should be observed here that equation (14) is obtained when treating the 
flow past an infinite yawed cylinder by the Pohlhausen method: for the external 
velocity components are U= ,(£) and V=constant and the application of the 
Pohlhausen method leads at once to equations (11) and (12). 


When p and k are known the velocities are found by expanding in powers 
of a:— 
u=u,+ au, +O (2?) 


v=2v,+O(a’) 
mx x0 { (ax- +x 


{xO 


=n 
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where y= { } 


In determining the skin friction it is observed that 


+0 (2°); 


_ 


The method has been applied to the flow considered in Part I, i.e. x(6)=€". In 
this case A, is independent of € when m is given and Cooke’ has shown that p is 
then also independent of £. He gives values of p for some values of 8 (=2m/(m+ 1)). 
When equation (15) is solved it is found that 


k CE —m+1 
where C depends only on £. 


The quantities 


—(3m—1) 1) 
= (m4+1) 


(3) =)... =c(2-%),/(E) =F" 


—(m—1) 


are comparable with F,” (0), F,” (0), G,” (0) of Part I and are given in Table Il 
together with F,” (0), F,” (0), G,” (0), A, and q(=1/p). When 8=0, A,=0 and 
the equations as given cannot be used; but the method is simply modified in this 
case and Table II includes values for 8=0. 


It will be observed that the best value of F,” (app.) occurs for a value of 8 
between 0-1 and 0-2 and, as is well known, the agreement between F,,” (app.) and 
F,” (0) is good for all positive 8. F,” (app.) has its best value at some slightly 
negative value of 6 and the error increases as 8 increases. At B=1 it is 16-5 per 
cent. of F,” (0). The agreement between G,” (app.) and G,” (0) is good for ~ 
positive and improves as £ increases. 


In each case as 8 becomes negative and tends to the value —0-1988, the error 
increases rapidly. 
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The other flow considered has x (€)=1-—€. In this case the Pohlhausen method 
applied to the corresponding two-dimensional flow is known to give bad results 
and so a good approximation is not to be expected from the present method, which 
depends on the accuracy of the two-dimensional approximation: but the figures 
given in Table III should indicate the trend of the three-dimensional effects. 

At the two-dimensional separation point, where A,= — 12 and £=0-156, p=0-949, 
k=0-25, 


vi 


vt 
55 


approximately while the ratios 


ME), 


have become infinite. 


It must be pointed out that the figures in Table III are not very accurate, 
especially for the larger values of €. More accurate values could be obtained only 
by using a smaller interval of € in solving the differential equations (14) and (15). 


It can be shown, by using a method similar to that used in Part I, that there 


is an exact solution with 
uy =f, (z2)+O(€) 


where z=4 (é~!»v~! and u, is the Howarth solution” for the two-dimensional flow 
U=1- 
From the expressions (16), (17) and (18) it can be shown that the Pohlhausen 


solution is of the form 
u,=F, (z,)+O(§) 


u,=4 F, (z,)+ 
v,=0{F, (z,)+O 
where z,=acé -1y-4 and a is a certain constant (—0-1714). 


and 


The same sbi are obtained hai the exact solution, so that it seems 
that near the origin of € the Pohlhausen method gives results in good agreement 
with the exact solution. It is impossible to determine at what value of € the 
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Pohlhausen and the exact solution begin to diverge seriously, but it is known that 
when £=0-120 (which corresponds to two-dimensional separation in the exact 
solution) the Pohlhausen solution for the two-dimensional flow gives poor results: 
and, with the three-dimensional solution depending on the two-dimensional, the 
Pohlhausen solution must be in serious error for the larger values of £. However, 
it seems probable that the Pohlhausen solution will indicate the general trend of 
the three-dimensional effects and it will be seen from the figures in Table III that 
the skin friction is everywhere increased so that, although no certain conclusion 
can be drawn about the effect on separation, it seems reasonable to suppose that it 
will be delayed by the three-dimensional motion. 


The figures do not indicate a large movement of the separation point as a 
consequence of the cross flow. To move the separation point from ¢=0-°156 to 
£=0-160 would require a value of 2 of about 0-12. 


Conclusion 


In discussing the three-dimensional boundary layer flow defined by the external 
irrotational velocity components U= y(£), V=ay it has been shown that in the 
particular cases for which x (¢)=é" and x(é)=1-€ the skin friction is increased 
as the result of the cross flow, and although the breakdown of the method of 
solution for the critical value of m in one case, and the inadequacy of the solution 
for large value of € in the other, make it possible to draw only tentative conclusions 
about separation, it is likely that this will be delayed. Although it does not seem 
possible te give a complete physical explanation of these results at the present stage 
it is clear that there is a transfer of energy from the weaker cross flow to the stronger 
flow in the main stream direction which enables the flow to continue farther before 
the retarding action of the boundary causes separation. 
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